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These lectures, given at the 2014 Theoretical Advanced Study Institute (TASI), 
are an introduction to what we know at present about dark matter and the ma¬ 
jor current experimental and observational efforts to identify what it consists of. 
They attempt to present the complexities of the subject, making clear common 
simplifying assumptions, to better understand the reach of dark matter searches. 


1. Introduction 

Observations on all astrophysical and cosmological scales, from the scale 
of dwarf galaxies up to the largest scales, indicate that dark matter (DM) 
constitutes about 25% of the content of the Universe. It cannot consist of 
atomic matter, which makes up stars, planets and ourselves. This known 
form of matter accounts for at most 5% of the content of the Universe. The 
remaining 70% consists of dark energy, a component which, unlike matter, 
has repulsive gravitational interactions. 

The nature of DM remains one of the most important open problems 
in science, and the efforts to resolve it are numerous. The hunt for DM is 
multi-pronged and interdisciplinary, involving cosmology and astrophysics, 
particle physics, direct and indirect detection experiments and searches at 
particle colliders. 

Many DM candidates have been proposed, many of them associated 
with beyond the Standard Model (SM) physics at the electroweak scale. 
Much has been said about the lightest supersymmetric particle, such as 
the lightest “neutralino”, or other Weakly Interacting Massive Particles 
(WIMPs) as “natural” DM particle candidates based on the “WIMP Mira¬ 
cle” argument (see e.g. [1]) This has lead to some degree of disillusionment 
about the prospect of discovering the nature of DM, based on the lack of 
any indication of new particle physics below energies of a few TeV at the 
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Large Hadron Collider (LHC), and up to much higher energy scales in rare 
flavor physics processes. It is therefore important to avoid at this point 
oversimplifying assumptions in our presentation of what we know about 
the DM and what not having found it so far means. 

We will start by reviewing the evidence for DM and what we actu¬ 
ally know about it. This is what defines the DM problem. As we discuss 
the properties the DM must have, we will mention DM candidates which 
accommodate these properties in various different ways. Then we will con¬ 
centrate on particle DM candidates, how to compute the relic abundance 
of WIMPs, and how they are searched for both in direct and indirect de¬ 
tection experiments and at the LHC. We will discuss along the way recent 
DM hints in direct and indirect searches and DM particle candidates pro¬ 
posed to account for them. We will then mention briefly axions and sterile 
neutrinos. 


2. Evidence For Dark Matter 

There is evidence for the existence of much more matter than what can be 
assigned to visible matter at all scales from dwarf galaxies to the largest 
cosmological scales (see e.g. [2] for a recent review). The excess is what is 
called dark matter (DM). 

Galaxies are the building blocks of the present Universe. They range in 
mass from 10 9 M 0 (M 0 is a solar mass) for dwarf galaxies to about lO 13 Af 0 . 
Our own barred spiral galaxy, the Milky Way, has a mass of about 10 12 M 0 . 
So far 25 satellite dwarf galaxies have been observed bound to the Milky 
Way, within a radius of several 100’s of kpc from the galactic center (lpc 
(parsec) = 3.26 light-years). The closest of them is the Sagittarius Dwarf 
Galaxy, at 70 kpc from the center. Galaxies come in groups, clusters and 
superclusters. The Milky Way is part of the Local Group, together with 
two other large galaxies, Andromeda and Triangulum, and many smaller 
satellite galaxies within a radius of several Mpc. The Local Group is in the 
outskirts of the nearest cluster, the Virgo Cluster, whose center is 10 Mpc 
from the Milky Way. At the level of superclusters (scales of 10’s of Mpc) 
the structure of the Universe is better described in terms of filaments, walls 
and voids. The dominant component of all these structures is DM. 

Evidence for the existence of DM was first discovered by Fritz Zwicky 
in the 1930’s. Measurements of the velocity dispersion of galaxies in the 
Coma Cluster led him to the conclusion that they could not be bound to 
the cluster by the gravitational attraction of the visible matter (stars, gas, 
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dust) alone [3]. The same argument applied to the gas bound to clusters, 
as well as weak gravitational lensing measurements of the total amount of 
mass in clusters, indicate that DM constitutes 6 times the mass in visible 
matter [4]. However, the seminal 1937 paper [3] of Zwicky had only about 
ten citations in the first forty years. 

The existence of DM was rediscovered in the 1970’s, this time at galactic 
scales, by Vera Rubin and others. The flat rotation curves of disk galaxies, 
i.e. the constant rotational speed v as function of the radius r beyond the 
visible disk, indicates that the mass continues to grow with the radius [6]. 
If all the mass of the galaxy M would be concentrated at the disk, an 
object of mass m orbiting the galaxy beyond the disk would experience 
a gravitational force F = GMmn/r 2 = mv 2 /r, and v = ( GM/r) 1 ^ 2 would 
decrease as r 1 / 2 . A constant v requires the mass M to grow with r. Most of 
the matter of galaxies resides in spheroidal “dark haloes”. This evidence, as 
well as weak [7] and strong [8] gravitational lensing measurements, indicate 
that galaxies have more than 4 times the amount of mass in stars, gas and 
dust. Dwarf galaxies are the most DM-dominated systems known. 

On cosmological scales, several observations are combined with Gen¬ 
eral Relativity to determine the composition of the Universe. The density 
fraction flj = p,;/p c , is the ratio of the energy density of each component 
to the “critical density” p c , the value of the total energy density needed 
for the Universe to be spatially flat, fi = has been measured to be 

very close to 11 = 1 [9]. The latest observations by the Planck collabo¬ 
ration of the anisotropies in the Cosmic Microwave Background (CMB) 
combined with others (see Table 10 of [9]) have lead to a composition of 
fl DM h 2 = 0.1187 ± 0.0017 of DM, n b h 2 = 0.02214 ± 0.00024 of ordinary 
matter (baryons, i.e. protons and neutrons) and Qde = 0.692 ± 0.010 of 
dark energy. Here, h = 0.697 ± 0.024 is the present expansion rate of the 
Universe, the Hubble constant, in units of 100 km/Mpc s (we quote 68% 
limits). Five independent measurements of the abundance of atomic matter 
in the Universe show that it amounts to less than 5% of the content of the 
Universe: the X-ray emission from galaxy clusters, the relative height of the 
odd and even peaks in the angular power spectrum of CMB anisotropies, 
the abundance of light chemical elements generated in Big Bang Nucle¬ 
osynthesis (BBN), Baryon Acoustic Oscillations and absorption lines of the 
light of Quasars. 

All measurements so far have confirmed the Big-Bang model of a hot 
early Universe expanding adiabatically for most of its lifetime of tjj = 13.798 
±0.037 x 10 9 y [9]. The CMB was emitted 379 ky after the Bang, when 
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the temperature of the Universe was T ~ 3 eV and atoms became stable 
for the first time. After this moment, called “recombination”, the Universe 
was populated by neutral constituents (the atoms) instead of plasma, and 
the mean free path of photons became very large on cosmological scales. 
Shortly before the emission of the CMB (when the Universe was about 100 
ky old and T ~ 1 eV) the Universe passed from being radiation-dominated 
to being matter-dominated. 

BBN is the earliest episode in the history of the Universe from which 
we have data, consisting of the relative abundance of the light elements 
produced then: D, 4 He and 7 Li. It took place between 3 and 20 minutes 
after the Bang, when the temperature of the Universe was T ~ MeV. In 
order for BBN and all the subsequent history of the Universe to proceed 
as we know it, it is enough that the earliest and highest temperature of 
the radiation-dominated period in which BBN happens is just larger than 
4 MeV [10]. 

Inflation is an early period of exponential expansion in the history of the 
Universe, advocated to produce through quantum fluctuations of a scalar 
field the small density inhomogeneities which seed the large scale structure 
of the Universe During inflation, the radiation and matter densities become 
zero and the Universe is repopulated during the “reheating” period, after 
inflation. The earliest (highest) temperature of the radiation-dominated 
period after reheating is called the “reheating temperature”, Trh ■ 

The recent BICEP2/Keck Array and Plank [11] upper limit on gravity 
waves produced during inflation (which replaced the BICEP2 earlier claim 
of a measurement [12]) determines an upper limit on Trh- The upper limit 
on the ratio of tensor gravitational waves and scalar density perturbations 
in the primordial plasma tells us that the potential energy density during 
inflation was V < (10 16 GeV) 4 . Under the assumption of an instantaneous 
reheating, the most efficient type of reheating, this energy goes completely 
into radiation, V — T^ H . If the reheating is not instantaneous, the energy 
density decreases with the expansion of the Universe and can be much 
smaller when transferred to radiation. Thus Trh 10 16 GeV. How small 
Trh can be depends on the details of the reheating period [13]. Thus, the 
lower limit Trh > 4 MeV [10] imposed by BBN holds. 

We do not know the thermal history of the Universe before its temper¬ 
ature was 4 MeV, and most DM particle candidates are produced during 
this period. 
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3. What do we know about Dark Matter? 

We usually say that DM is neutral and stable, with very weak interactions. 
Let us review what we actually have learned about DM in the 80 years 
since we first found evidence of its existence. 

-1- DM has attractive gravitational interactions and is either sta¬ 
ble or has a lifetime tu- This is clear: DM is still present in the 
Universe and behaves as does regular matter for gravitational interactions. 

We have, in fact, no evidence that DM has any other interaction but 
gravity. Thus, one can wonder if the many observational evidences for 
DM are instead showing departures from the law of gravity itself. This 
general idea cannot be tested unless expressed in a particular model and the 
most successful model of this type proposed so far has been the “Modified 
Newtonian Dynamics” (MOND) [14], and its covariant version, the “Tensor- 
Vector-Scalar” (TeVeS) gravity model [15]. But they cannot replace DM, 
as described below. 

-2- MOND with only visible matter is not enough at scales larger 
than galactic scales. In 1983 M. Milgrom [14] proposed MOND as an 
alternative explanation for the flat rotation curves of galaxies. In MOND 
F = ma becomes F = g{a)ma, where g(a) deviates from unity only for 
very small accelerations a<«o, for which g = a/ao- Thus the gravitational 
force acting on a body of mass in in an orbit of radius r around a galaxy 
of mass M leads to F = GMm/r 2 = mag, which for large r, i.e. a <C do, 
yields a = \JGMa 0 /r. Equating this with the centripetal acceleration 
a = v 2 /r, one gets a constant orbital speed v = (GMao) 1 / 4 . 

MOND is in good agreement with galaxy-scale observations for ao — 
1.2 x 1(T 10 m/s 2 , without any need for DM [16] (and for this a o, the effects 
of MOND are too small to be measurable in laboratory or solar-system scale 
experiments). But it fails at larger scales [16], in particular in the “Bullet 
Cluster” [17], unless some form of DM is introduced [18]. MOND is only 
a non-relativistic theory, thus it cannot be used where General Relativity 
is needed. TeVeS [15] is a relativistically covariant theory in which the 
tensor field in Einstein’s theory of gravity is replaced by scalar, vector and 
tensor fields, which interact in such a way to yield MOND in the weak-field 
non-relativistic limit. 

In the “Bullet Cluster” [17], discovered in 2004, the baryonic matter 
(hot gas observed in X-rays) is at the center, spatially segregated from the 
two lateral gravitational potential wells (measured via weak gravitational 
lensing). The explanation of this system based on DM is that two galaxies 
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collided, leaving behind the interacting gas, while the DM of both galaxies 
passed through (which yields an upper limit on the self interaction of the 
DM). To explain this system TeVeS requires 2-3 times more matter than 
accounted for by the visible matter alone [16, 18]. TeVeS advocates propose 
the existence of some “cluster baryonic dark matter” (CBDM) [19], for 
example. But once some kind of extra matter is necessary, it then becomes 
an issue of taste which type of DM one considers more likely. 

-3- DM is not observed to interact with light. It is not observed to 
emit, reflect or absorb light of any frequency. This means that most of the 
DM must have a small enough electromagnetic coupling or be very heavy. 
Upper limits are derived from background light at all frequencies [20]. 

The DM could be neutral, maybe with a small electric or magnetic 
dipole moment [21, 22, 23] or an anapole moment [21, 24]. It could have 
a small effective electric charge, such as that of “Milli-Clrarged DM” [25]. 
Charged DM particles can also act as nearly collisionless if their mass is 
sufficiently large [26, 27]). 

“Millicharged DM” (see e.g. [25, 28]) can be part of a “dark sector” 
which couples to SM particles only though the admixture of a “dark pho¬ 
ton” of a “dark U’(l)” gauge symmetry with the usual photon (or the Z 
boson). This could occur through a “kinetic mixing” [29, 30] of the form 
eF llv F' liv , where F^ u and F are respectively the field-strength tensors 
for the hypercharge and dark U(l) gauge groups. After diagonalizing the 
kinetic terms, the usual photon has an admixture e of the dark photon. 
Thus if the DM particle couples with charge Q' to the dark photon, the 
effective coupling with the usual photons is eQ' and several limits imply 
that e < 1CT 3 [28, 31], 

“Millicharged DM” appears naturally as part of a complicated dark 
sector. It could be “Atomic DM”, with dark protons and dark electrons 
forming dark atoms [32, 33, 34, 28], or “Mirror DM” [35] whose Lagrangian 
is a copy of that of the SM, but for the mirror particles (see e.g. [45] 
and references therein). These are some possibilities in which the dark or 
“secluded” sector has some of the richness of the visible sector, with hidden 
gauge interactions and flavor (e.g. [36] and [37]). 

Observational upper limits on the cross section of elastic DM -photon 
interactions can be surprisingly large, e.g. a recent limit from simulations 
of Milky Way sub-haloes is cr|]|| t i c 7 < 4 x 1CT 33 cm 2 (rn/GeV) [38], 

An important consequence of the small interaction of DM with light is 
that the DM cannot cool by radiating photons during galaxy formation. 
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-4- The bulk of the DM must be dissipationless, but part of it 
could be dissipative. “Dissipationless” mean that the DM cannot cool by 
radiating as baryons do to collapse in the center of disk galaxies. Otherwise, 
their extended dark halos would not exist. Galaxies start as structures made 
of the primordial admixture of dark and visible matter. Then, visible matter 
dissipates energy by emitting photons and falls into the potential well of 
the object. Because this emission is isotropic, the visible matter preserves 
any angular momentum it might initially have. Thus as it collapses to the 
center, it increases its angular speed until it becomes unstable towards the 
formation of a disk, which thus rotates much faster than the dark halo. 

While most of the DM must be nearly dissipationless, a small fraction of 
it could be dissipative. Part of the matter in the galaxy, the visible matter, 
is dissipative and its presence does not disrupt the stability of dark haloes. 
Thus a similar fraction, 5-10%, of the DM could also be dissipative and 
even form a “Dark Disk”. This is the idea behind “Partially Interacting 
DM” (PIDM) and “Double Disk DM” (DDDM) [39]. The dissipative DM 
component could emit “dark photons” or other “dark” particles. 

A Dark Disk was shown to arise in some simulations of galaxy formation 
including baryonic matter besides the usual non-dissipative Cold DM [40], 
but with dissipative DM it should be a pervasive feature of all disk galaxies. 
-5- The mass m of the major component of the DM has only been 
constrained within some 80 orders of magnitude. There is a firm 
upper bound m < 2x 10 -9 M 0 = 2x 10 48 GeV at the 95% CL. It comes from 
unsuccessful searches for MACHOS (“Massive Astrophysical Compact Halo 
Objects”) in the dark halo of our galaxy using gravitational microlensing 
with the Kepler satellite [41], and from the ground-based MACHO and 
EROS surveys [42], combined with bounds on the granularity of the DM 
for masses larger than 30 M 0 [43]. Microlensing is a type of gravitational 
lensing in which the multiple images of the lensed star are superposed, 
producing a magnification of the star flux if an object passes near the line 
of sight to the star as it moves through the dark halo. Above this limit, 
MACHOS can account for only a small fraction of the dark halo of our 
galaxy. 

Among the best candidates for MACHOS are “Primordial Black Holes” 
(PBH) [44]. These are hypothetical black holes that could be created in a 
primordial phase transition, maybe during inflation (see e.g. [45, 46]). How¬ 
ever, several limits apply to PBH which do not apply to other MACHOS and 
constrain the fraction of DM in PBH to be < 1 for almost any TOpbh mass. 
PBH lighter than 10 15 g would have decayed by Hawking radiation. Not 
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having observed this radiation in 7 -rays excludes rapHB<10 17 g. Accretion 
of PBH during star formation which would destroy neutron stars excludes 
10 16 g<mpHB<10 22 g. Part of this last range, 10 17 g<mpHB<10 2 °g, is also 
excluded by the non-observation of “femtolensing” (in which the lensed 
images are separated by femto arc-seconds) of Gamma-Ray Bursts. The 
only narrow window remaining for PBH to make up all the DM, 10 22 g to 
4 x 10 24 g (i.e. 6 x 10 45 GeV to 2 x 10 48 GeV), is being challenged by less- 
certain excluding arguments. For a current review on PBH and references 
to all these limits see [46]. PHB DM does not require physics beyond the 
SM (except some special inflation potentials). 

There could be other DM candidates which might arise within the SM, 
generically called “Macroscopic DM” [27]. How to generate enough of 
them to account for all the observed DM is problematic. These include 
“Nuclearites”, “Strangelets”, “Strange Baryon Q-Balls”, “Baryonic Color 
Superconductors”, “Compact Composite Objects”, “Strange Chiral Liq¬ 
uids Drops”, and “Compact Ultradense Objects”. Macroscopic candidates 
could arise in physics beyond the SM too, e.g. “Supersymmetric Q-balls”. 
Various limits on them are presented in [27] (see references therein for all 
these candidates). 

The limits just presented, and the fact that particle candidates can have 
the right relic abundance to be the DM, constitute the only observational 
arguments we have in favor of elementary-particle candidates for DM. 

The lower limit on the DM particle mass is less well determined. It is at 
least 10 -31 GeV since there is a concrete particle candidate proposed with 
this mass, “Fuzzy DM”. This is a boson with a de Broglie wavelength of 1 
kpc [47]. For particles that reached thermal equilibrium in the early Uni¬ 
verse the lower limit is 0.2-0.7 keV instead, due the maximum occupation 
number these particles can have in different structures in the Universe [48]. 
-6- DM has been mostly assumed to be collisionless, however the 
upper limit on DM self-interactions is very large. The best limits 
on self-interactions are derived from the “Bullet Cluster” [17] (where the 
DM of the two colliding galaxies had to pass through each other) and the 
non-sphericity of the halos of galaxies and galaxy clusters. The upper 
limit is huge: a se \i/m < lcm 2 /g ~ 2 barn/GeV ~ 2xl0 -24 cm 2 /GeV. By 
comparison, the neutron capture cross section of uranium is a few barns! 
The limit on the cr S eif /m ratio comes from requiring that the self-interaction 
mean free path, X m f p ~ l/n< 7 se if = m/pa se if, be long enough (n = p/m is 
the DM number density, and the density p of the system is measured). 



Explaining the existence of DM cores in dwarf galaxies may constitute 
a problem for the usual collisionless Cold DM models. Collisions of the DM 
particles with themselves would erase small scale structure and turn cuspy 
central density profiles into cored profiles in dwarf galaxies, if the cross 
section is close to the mentioned upper limit [49]. DM with this large self¬ 
interaction is called “Self-Interacting DM” (SIDM) [50]. DM with less than 
an order of magnitude smaller <r se if is indistinguishable from collissionless. 

With SIDM the reduced central densities of dwarf galaxies imply re¬ 
duced velocity dispersions, which can alleviate [51] the “too big to fail 
problem” [52] of collisionless CDM (see below). There are many particle 
models for SIDM, most with complicated “dark sectors” (see e.g. [45]) 

-7- The hulk of the DM is Cold or Warm, thus particle DM re¬ 
quires physics beyond the SM. The DM is classified as hot, cold or 
warm according to how relativistic it is when galactic-size perturbations 
enter into the horizon (i.e. when these perturbations become encompassed 
by the growing horizon ~ ct). This happens when T — keV. Hot DM 
(HDM) is relativistic, Cold DM (CMD) is non-relativistic and Warm DM 
(WDM) is becoming non-relativistic at this moment (see e.g. [53]). 

The presence of DM in the early Universe is necessary for the formation 
of structure in the Universe. Structures in baryons cannot grow until re¬ 
combination, when atoms become stable. Before then the photon pressure 
in the plasma prevents it. But at recombination baryons must fall into al¬ 
ready formed potential wells of DM, or there would not be enough time to 
form the structures we observe now. Thus, it is the DM which determines 
the major features of the large-scale structure of the Universe. 

Perturbations in the DM survive at horizon crossing only if the DM 
is non-relativistic. With HDM, primordial galaxy-size density inhomo¬ 
geneities would not survive, superclusters would form first and later galax¬ 
ies through fragmentation. This does not reproduce the observed Universe. 
With CDM, inhomogeneities much smaller than galaxy size survive, thus 
galaxies and clusters incorporate many smaller structures which form first. 
Some of them are not entirely tidally destroyed, thus, a large number of 
substructure is expected within CDM dark haloes. With WDM the smaller 
structures formed first are of the size of dwarf galaxy cores, thus there is 
much less substructure within large haloes than with CDM (see e.g. [53]). 

Either CDM or WDM can account for all the large scale structure ob¬ 
servations. The difference between them is at the dwarf-galaxy scale, where 
observations and their interpretation are still not conclusive. 

Very high resolution simulations of structure formation in the Universe 
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with only CDM (assumed to be collisionless) find within a halo similar 
to that of the Milky Way, of the order of 10 subhaloes so massive and 
dense that they seem “too big to fail” to form lots of stars within them 
and be visible. But the Milky Way and Andromeda do not seem to have 
satellite galaxies with stars moving as fast as would be expected in these 
dense sub-haloes. This constitute the so called “too big to fail” problem of 
CDM [52]. Due to the large uncertainty in the mass of the Milky Way, this 
problem may be alleviated in our galaxy if its mass is sufficiently small, 
but the problem persists in Andromeda and in the Local Group [54]. There 
are indications of the existence of a core (a region of constant density) in 
the DM density profile of dwarf galaxies, instead of the cusp (with density 
growing towards the center) predicted by CDM-only simulations. 

Both these potential problems of (collissionless) CDM may disappear 
once the effect of visible matter is fully taken into account, or with WDM 
or SIDM (Warm or Self-Interacting DM) [50, 49, 51] instead of CDM (see 
also [45] and references therein). With WDM too few satellite galaxies 
might be produced in the Milky Way, if its mass is small enough [55]. 

The only particles in the SM which are part of the DM are neutrinos. 
They are lighter than 1 eV and remain in thermal equilibrium in the early 
Universe until temperatures of a few MeV (see below), thus they are HDM. 
There are no CDM or WDM particle candidates in the SM but there are 
many in extensions of the SM. Candidates for CDM are “axions”, WIMPs, 
and many others. Sterile neutrinos and some non-thermal WIMPs (see 
below), among others, are good WDM candidates. 

Because of spontaneous symmetry breaking arguments, completely in¬ 
dependently of the DM issue, we do expect new physics beyond the SM to 
appear at the electroweak scale. Many extensions of the SM were proposed 
because of this reason (supersymmetry, technicolor, large extra spatial di¬ 
mensions, the “Little Higgs” model, the Inert Doublet model etc.), which 
provide both the main potential discoveries at the LHC, and also DM can¬ 
didates. However, the DM physics may be entirely different from that at 
the electroweak scale. Many models have been proposed in recent years just 
to account for DM hints in direct and indirect DM searches. We already 
mentioned some of them, and will mention more later in the lectures. 

When computing the properties of DM candidates we should keep in 
mind our assumptions, among them our cosmological assumptions. 

-8- Most DM candidates are relics from the pre-BBN era, from 
which we have no data. The computation of the relic abundance and 
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primordial velocity distribution of particle DM candidates produced before 
the BBN temperature limit of 4 MeV depends on assumptions made regard¬ 
ing the thermal history of the Universe. With different viable cosmological 
assumptions (see e.g. [56]), the relic density and velocity distribution of the 
DM candidates may change considerably (see section 4.1.4). 

As already pointed out early on (e.g. by J. Barrow in 1982 [57]), if 
we ever discover a DM particle from the pre-BBN era, we would use its 
properties as a cosmological probe to gain information about this epoch. 


4. The relic abundance of Dark Matter particle candidates 

We usually characterize DM particle candidates according to how they are 
produced as “thermal” or “non-thermal” relics (for a review see e.g. [56]). 
“Thermal” relics are produced via interactions with the thermal bath, reach 
equilibrium with the bath and then “decouple” or “freeze-out” when their 
interactions cannot keep up with the expansion of the Universe. Chemical 
equilibrium is achieved when reactions that change the number of parti¬ 
cles are faster than the expansion rate of the Universe H (or the reaction 
time is shorter than the lifetime of the Universe ~ H~ 1 ). After chemical 
decoupling or freeze-out the number of particles per comoving volume re¬ 
mains constant. After kinetic decoupling, the exchange of momentum with 
the radiation bath ceases to be effective. “Non-thermal” DM particles are 
all those not produced in this way. For example, they could be produced 
via the decay of other particles, which themselves may or may not have a 
thermal abundance. 


4.1. Thermal relics 

The evolution of the number density n of thermal DM particles x and 
antiparticles x interacting with particles in a thermal bath is given by the 
Bolzmann Transport Equation (see e.g. [53]), 

^ =-3Hn-(a A v) (n 2 (1) 

assuming for simplicity that x and X can only annihilate and be created in 
XX pairs and that there is no particle-antiparticle asymmetry, n x = n x = n. 
Both of these assumptions arc fulfilled for the neutralino and many other 
candidates that coincide with their own antiparticles. Here t is time, H 
is the expansion rate of the Universe (the Hubble parameter, H = a/a, 
where a is the scale factor of the Universe), a a is the XX annihilation 
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cross section, and v is the relative speed of the annihilating DM particles. 
The angle brackets ( ) denote an average over the thermal distribution of 
momenta of the DM particles (an average over initial states and sum over 
all final states, see e.g. [58]), and tieq is the equilibrium number density at 
the temperature T of the thermal bath. 

The three terms on the right-hand side of Eq. (1) account respectively 
for the expansion of the Universe, particle annihilation and creation (see 
e.g. [53]). If particles were neither created nor destroyed, the solution n ~ 
a~ 3 would only be diluted by the increase of a, with dn/dtt = —3 (a/a)n = 
—3Hn. With only annihilation, n ~ e~ t ^ iA and dn/dt = — n/tA , where 
the typical annihilation time f/t is proportional to the mean free path t>A — 
A mfp/v = l/(tr J 4 ?rn). Thus the term proportional to n 2 accounts for the 
annihilation. Include the creation, stop the expansion when the thermal 
bath temperature is T and wait until equilibrium is achieved, in which case 
n = 71eq and dn/dt = 0 , thus the creation term must cancel the annihilation 
term with n = Ueq. 

Consider the entropy density s = (27r 2 /45 )g s - e s(T) T 3 , where g s - eS(T) 
is the effective entropy number of degrees of freedom [53, 59], and T the 
photon temperature. Eq. (1) and the conservation of entropy per comoving 
volume, S = sa 3 =constant, which implies ds/dt = — 3Hs , can be combined 
into a single equation for the dimensionless variables Y = n/s and x = 
m/T, 


( 2 ) 


dY Ids 2 \r2 \ 

~dx-Mdx {(TAV){Y _Yb « } - 
Notice that with S = sa 3 constant, Y is proportional to the number of par¬ 
ticles per comoving volume na 3 . When g s - e s{T) is approximately constant 
(thus a ~ T -1 ), a good approximation most of the time, Eq. (2) becomes 


J^dY = _TA \ ' 

Yeq dx H 1 [y£q) 


( 3 ) 


where = tieq (&v) is the equilibrium annihilation rate. Eq. (3) clearly 
shows that the number of particles per comoving volume becomes constant 
(dY/dx = 0) when Ta/H « 1. At high T, Ta > H for particles that 
are in equilibrium, but T^ decreases with decreasing T faster than H , and 
crosses H at the chemical decoupling or freeze-out temperature T = Tf a , 


Ta(T} 0 ) — ( <J Av) T _ Tfo nEQ(Tfo) — H(Tf 0 ). (4) 

H is related to the total energy density of the Universe p by Friedmann’s 
equation H = \J8irGp/3 [53, 59]. G is the gravitational constant and 


12 



G = Mp 2 defines the Planck mass Mp ~ 10 19 GeV. For a radiation- 
dominated Universe, p ~ T 4 and thus H ~ T 2 /Mp. 

Let us estimate X/ 0 for particles that are relativistic (e.g. active neutri¬ 
nos) or non-relativistic (e.g. WIMPs) at decoupling, using estimates of the 
annihilation cross section into light SM fermions of mass rrif « T via the 
exchange of a mediator of mass M and coupling g. 

For relativistic particles, in < T, thus T is the dominant quantity of 
mass dimension which can appear in the numerator of the cross section to 
compensate the 1/M 4 propagator factor (the mass dimension of a is —2) 


n 4 

R y t>2 

aA -W T - 


(5) 


For weak interactions g 4 /M 4 ~ G 2 F ( Gf ~ 10 5 /GeV 2 is the Fermi con¬ 
stant). For non-relativistic DM particles there are two possible situations. 
If M > m > T, m and not T appears in the cross section, namely 


nr ^ 9 2 

A " M 4 ™ ■ 


( 6 ) 


If instead the DM mass is much larger than any other quantity with mass 
dimension, including M, only m appears in the cross section, 


NR 

a A 



(7) 


4.1.1. Particles relativistic at decoupling - active neutrinos 
Active neutrinos interact via weak interactions, and Eq. (4) becomes 

rp 2 

T A (T fo ) cs n EQ (T fo )a%(T fo ) ~ G%T% ~ H(T fo ) ~ (8) 

We used that for relativistic particles u-eq — T 3 , and that Eq. (5) gives 
a)/ ~ GpT 2 . We have also used the form of H valid for a radiation- 
dominated Universe. Substituting in the numerical values for Gf and Mp 
we get Tf 0 ~ few MeV, the same result as from a full calculation. Since the 
active neutrinos have masses smaller than a few eV, they are relativistic at 
decoupling. These thermal neutrinos constitute the Cosmic Neutrino Back¬ 
ground and are the most abundant particles in the Universe after the CMB 
photons. Their temperature is smaller than the CMB photon temperature, 
T„ = (4/H) 1 / 3 ^, because e+e~ pairs annihilate after neutrinos decouple, 
so that their entropy goes almost entirely to photons. 
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4.1.2. Particles non-relativistic at decoupling- the “WIMP Miracle” 

For particles with m » Tf a , ueq = g x ( m T /27t) 3 ^ 2 e~ m ! T where g x is the 
X number of degrees of freedom (normalized so that g^ = 2 for a photon). 
Thus Eq. (4) becomes 


r(T /o ) cs aT n E Q ^ (T fo 


( mT fc 

I 2tt 


3/2 J-i2 

e -m/T So „ H{Tfo) „ 


(9) 

The crossing of the terms proportional to e~ m l T i° and T 2 a depends mostly 
on the exponential factor: it will happen when the argument of the expo¬ 
nential is neither very small nor very large, i.e. m/Tf 0 ~ 1 (actually one 
finds m/Tf 0 of 0(10)). Using Tf a ~ m in Eq. (9) we get 


«EQ (Tfo) ~ 


T 2 
1 So 


rNR/ 


( 10 ) 


'A l (T fo ) aT(T fo ) 

After decoupling, at T < Tf 0 , the number density n(T) only decreases as 
T 3 due to the expansion of the Universe, n(T) = nEQ(7/ 0 )(T 3 /T 3 0 ). Thus, 
using Eq. (10), the DM density at T < Tf a is 

J~<3 


p{T) = to n(T) ~ 


rNR/ 


( 11 ) 


’ A \T fo )m? aT{T fo y 
With this simple argument [60], we obtain, the crucial result that the relic 
density is inversely proportional to the annihilation cross section at Tf 0 . 

The numerical solution of the Boltzmann equation, assuming the Uni¬ 
verse is radiation dominated and there is conservation of entropy in matter 
and radiation, is illustrated in Fig. 1. At high temperatures Y closely tracks 
its equilibrium value until at freeze-out, when annihilation and creation of 
WIMPs become ineffective. Then, Y becomes constant. The freeze-out for 
weak-strength interactions occurs at Xf 0 = m/Tf 0 ~ 20, when the typical 
WIMP speed is Vf a = ( 2>Tf 0 /2m) 1 / 2 ~ 0.27c, and the relic density is 


flh 2 ~ 0.1 


Xf 0 \ f 60 3 x 10 26 cm 3 /s 


(t)( 


9eft ) 


36/ Xf 0 


( 12 ) 


where (cr^ R c) ~ a + b(v 2 ) + 0(v 4 ) (a and bv 2 correspond to s-wave and 
p-wave annihilation, respectively). g e g is the effective number of degrees 
of freedom in the radiation bath at freeze-out (which coincides with the 
number g s - e g in the entropy density except at very low temperatures). 
Within the SM, g eS is 10.75 at IMeV < T < 100 MeV, is about 60 above 
the QCD phase transition and up to T ~ GeV, becomes about 90 at T ~100 
GeV, and reaches 106 above the top quark mass [53, 58, 59]. 
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Figure 1. Typical evolution of the WIMP number density per comoving volume Y in 
the early Universe assuming no particle-antiparticle asymmetry as function of x = m/T. 
Here Y eq is the equilibrium value of Y. Fig. from [56]. 



lg [eVj 


Figure 2. Relic density €l v h? of a 4th generation thermal Dirac neutrino with SM 
interactions as a function of the neutrino mass m v (solid line). Fig. from [56]. 


Smaller annihilation cross sections lead to larger relic densities (“the 
weaker prevail”). WIMPs with stronger interactions remain in equilibrium 
longer, hence decouple when their density is further suppressed by a smaller 
Boltzmann factor (e.g. for strong interactions, Xf Q ~ 45 instead of 20). 

The “WIMP miracle” refers to the fact that for typical cross sections 
of weak order of magnitude, Eq. (12) gives the right order of magnitude of 
the DM density f Iti 2 ~ 0.1 for WIMP masses in the GeVs-to-TeV range. 
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This is shown in Fig. 2 for the original WIMP, studied in 1977 by Lee 
and Weinberg [61], an active Dirac neutrino of a 4th generation with no 
lepton asymmetry. Although this candidate has been rejected by a combi¬ 
nation of LEP and direct detection limits, it is still interesting as a ped¬ 
agogical example. Fig. 2 shows its relic density fI„/i 2 as function of its 
mass m „. For to„ < MeV these neutrinos are relativistic at decoupling, 
thus n v is the same for any mass and f l„h 2 ~ m„. For m „ > MeV they 
are non-relativistic at decoupling and Eq. (12) applies. Their annihilation 
is mediated by the Z-boson, thus the Z-resonance in the cross section at 
m v ~ Mz /2 gives rise to the characteristic V shape in the fl„/i 2 curve. For 
m v < M z , ( UannihV) — G 2 F ml ~ {m v / GeV) 2 1CT 26 cm 3 /s (Eq. (6)) and 
for m v > M z , {<J annih v) ~ a/m 2 ~ (m„/TeV) -2 10~ 26 cm 3 /s (Eq. (7)) . 
Here, a is the electromagnetic coupling constant. This shows that GeV- 
and TeV-mass WIMPs have the relic abundance of the DM 

Above m„ ~ 100 GeV, new annihilation channels into Z or W boson 
pairs open up, and when they become dominant, the annihilation cross 
section increases. This decreases D„/i 2 . For larger ?n„, a general unitarity 
argument [62] limits f \ v h 2 to be below the dashed line in Fig. 2, and then the 
neutrino Yukawa coupling becomes too large for perturbative calculations 
to be reliable (indicated by a question mark in Fig. 2). 

The relic density of other WIMP candidates exhibits features similar 
to that of the Dirac neutrino just discussed. State-of-the-art calculations 
of WIMP relic densities strive to achieve a precision comparable to that of 
observations, around a few percent, given the assumptions made. Notice 
the assumptions that go into the calculations presented above. 1) The DM 
is in thermal equilibrium and then decouples while the Universe is radia¬ 
tion dominated. 2) There is no entropy change in radiation plus matter, 
either during or after decoupling. 3) There is no DM particle-antiparticle 
asymmetry. 4) The DM particles are stable. A change in any of these as¬ 
sumptions could lead to a very different relic abundance: the DM could be 
asymmetric, the Universe may not be radiation dominated at decoupling, 
the entropy per comoving volume of radiation plus matter could change, 
DM WIMPs could be produced in the decay of other particles, thermal 
or not, or WIMPs could be unstable and decay into the DM after they 
decouple. 

This last possibility is the “SuperWIMP” scenario [63], in which the 
DM particles could have interactions much weaker than weak (super-weak) 
because they inherit the correct DM density from unstable thermal WIMPs 
which decay to them at late times. 
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4.1.3. Asymmetric Dark Matter (ADM) 

We owe our very existence to a particle-antiparticle asymmetry in baryons 
so, why not also the DM? Baryons B (protons and neutrons) decouple 
while they are non-relativistic. If the number of B and B would be the 
same before and during their decoupling, then after decoupling ns/s = 
Ub/s ~ 1(T 19 , and thus f Ib = Dp ~ 10 -11 (see e.g. [53]). What we observe 
instead is Db — 0.05 and = 0- Thus an early Baryon Asymmetry 
Ab = {tib — n§)/s — 1CT 10 must exist, so that ns cannot be smaller than 
sAb, as ng goes to 0. 

The idea of ADM is almost as old as the “WIMP Miracle”. It was 
proposed first in 1985 [64] for candidates with mass of 100’s of GeV in the 
context of technicolor models, and in 1986 [65] for particles of a few GeV of 
mass (in the context of models for “Cosmions”, DM particles which could 
cool the central region of the Sun). The idea has been pursued ever since, 
and became very popular for 1-10 GeV mass WIMPs in recent years (see 
e.g. [66] and references therein). The simplest models [64, 65] assume a 
similar asymmetry of the DM and baryons, Abm — Ab, thus their relic 
number densities are similar, udm — tib and 

ftpM ^ n DM m ^ ^ m 

Db tibitib tub lGeV 


If m ~ 5 GeV we obtain the right ratio SIdm/^b — 5. Thus ADM explains 
the present Ddm/^b ratio, which is otherwise a coincidence. 

Fig. 3 shows the evolution of Y x = n x /s and Y x = n x /s while particles 
X and antiparticles \ are in thermal equilibrium in the presence of an 
asymmetry Y — Y = A > 0. The equilibrium number densities of \ and X 
in this case differ by the chemical potential /i x (in equilibrium /i x = — fj, x ): 


n EQ _ 


= 9x 




EQ 

L X 


(y“) 7 e (_m_ ' 1 * )/T . (14) 


Therefore n x ®(n x = 0) [exp (n x /T) — exp (— jj, x /T)\ = As. Notice that 
V X E Q cannot become smaller than A (it approaches A as Y^ goes to zero). 
V E< ^ decreases exponentially until \ freezes-out, when ~ H , where 
~ n E< ^ ~ AT 3 . The previous Boltzmann Eqs. (1) and (2) become (if 
the DM particles annihilate and are created only in XX pairs) 

+ % Hn x = - ( a xx v ) ( n x n x ~ n x Qn x Q )• ( 15 ) 
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Figure 3. Typical evolution of the equilibrium number density per comoving volume 
in the early Universe of particles x an d antiparticles x, Y^^ and Y'?^, in the presence 
of an asymmetry Y — Y = A > 0, compared with Y E Q for no asymmetry (chemical 
potential /r = 0), as function of x = m/T. Fig. from [69]. 


and 


x dY x _ I\ ( Y x Yx 
Y^ Q dx H \ y^ Q Y eq 

X \ X x 


(16) 


where r x = (cr xx v} is the annihilation rate of particles X- The corre¬ 
sponding equations for \ are obtained by exchanging \ and X- Annihila¬ 
tions cease once the minority component (here \) decouple, when T x < H. 

In the standard cosmology (i.e. assuming radiation domination at de¬ 
coupling and entropy conservation) the relic density of the minority compo¬ 
nent [67, 68] is exponentially small with respect to the majority component 
density, which means that there is no DM annihilation after decoupling. 
This is one of the main characteristics ADM is assumed to have. However, 
this is a pre-BBN cosmology dependent feature. In some non-standard cos¬ 
mologies, the present ADM annihilation rate can be very large (see e.g. [69]). 


4.1.4. Non-Standard pre-BBN cosmologies 

The relic density and velocity distribution of many WIMPs and other DM 
candidates, e.g thermal WIMPs of mass m > 80 MeV, heavy sterile neu¬ 
trinos and axions, depend on the unknown characteristics of the Universe 
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before T ~ 4 MeV [10], when the content and expansion history of the 
Universe may differ from the standard assumptions. In non-standard pre- 
BBN cosmological models, the WIMP relic abundance may be higher or 
lower than the standard abundance. The density may be decreased by re¬ 
ducing the rate of thermal production (through a low T RH < Tf 0 ) or by 
producing radiation after freeze-out (entropy dilution). The density may 
also be increased by creating WIMPs from decays of particles or extended 
objects (non-thermal production) or by increasing the expansion rate of the 
Universe at the time of freeze-out (see e.g. [56, 70]). 

Not only the relic density of WIMPs but their relic velocity distribution 
can differ from the standard predictions. WIMPs could be warmer and even 
WDM [71], if they are produced in the decay of a heavy particle and do not 
exchange momentum with the thermal bath after, or could be colder [72]. 

Non-standard cosmological scenarios are more complicated than the 
standard scenario (e.g. to generate the baryon number asymmetry). They 
contain additional parameters that can be adjusted to modify the DM relic 
density. However these are due to physics at a high energy scale, and do 
not change the particle physics model at the electroweak scale or below. 


5. WIMP DM searches 

WIMP’s are actively searched for in direct and indirect DM detection ex¬ 
periments, and at colliders. 

Direct searches look for energy deposited within a detector by the DM 
particles in the dark halo of the Milky Way. They are sensitive to even 
a very subdominant WIMP component of the DM, if the scattering cross 
section as is large enough. For thermal WIMPs (p ~ 1 /a a), the event rate 
R ~ pas ~ as/a a- Because of the crossing symmetry relating annihilation 
and scattering amplitudes, if the annihilation cross section a a is large so is 
as- Therefore the rate may remain large even for WIMPs with a very small 
relic density (see e.g. [73]). Direct searches are subject to uncertainties in 
the local dark halo characteristics, and are relatively insensitive to DM that 
couples to leptons only (“leptophillic”). They would fail if the WIMP is so 
light its signal is below threshold, or if as is too small. 

Indirect searches look for annihilation (or decay) products of WIMPs 
in the dark haloes of the Milky Way and other galaxies, as well as WIMPs 
accumulated within the Sun or Earth. These searches are sensitive to inter¬ 
actions with all SM particles and directly probe the annihilation process in 
the “WIMP miracle”. The caveat to this type of search is that the DM may 
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not annihilate (e.g. if it is asymmetric and consists exclusively of particles 
and no antiparticles) or decay. Experimental sensitivities to several types 
of signal are expected to improve greatly in the coming decade, but the 
discovery of the DM through these searches requires understanding astro- 
physical backgrounds and the expected signal is subject to uncertainties in 
dark halo characteristics. 

At colliders, in particular the LHC, WIMPs are searched for as missing 
transverse energy. The caveats to this type of search is that the DM particle 
may be too heavy to be produced (above a few TeV at the LHC) or its 
signal may be hidden by backgrounds. Also, a signal produced by a particle 
escaping the detectors with lifetime ~ 100 ns cannot be distinguished from 
one with lifetime > 10 17 s as required for DM particles. Hadron colliders 
are relatively insensitive to DM that interacts only with leptons. 

All three types of searches are independent and complementary to each 
other. Even if the LHC finds a good DM candidate, to prove that it is 
the DM (and that it did not decay on cosmological time scales or was not 
produced in large enough amounts in the early Universe) we will need to 
find it in the dark haloes of our galaxy and other galaxies. 


5.1. Direct Detection of WIMPs 

The flux of WIMPs passing through a detector on Earth can be large, 
nv = pDMv/m ~ 10 7 (GeV/m)/cnr 2 s. However the expected energies de¬ 
posited and interaction rates are very small, keV to 10’s of keV and less 
than an event per 100 kg per day. Measuring these small energies and 
rates requires a constant fight against backgrounds. The experiments need 
to be underground, in mines or tunnels under mountains, to shield them 
from cosmic rays. The signal should consist of single hits and be uniform 
throughout the volume of the detector (this is a powerful way of discrimi¬ 
nating agains neutrons, which tend to interact multiple times in a detector 
and closer to its surface). 

Most direct searches are non-directional, but some, still in the stage 
of development, attempt to measure the recoil direction beside the energy 
(see e.g. [74]). In a directional search, 10’s of events would be enough to 
determine the direction of arrival and thus the DM origin of the signal. 
Without directionality, the unmistakable signature of DM is a few percent 
annual modulation of the rate, due to the variation in the velocity of the 
“WIMP-wind” on Earth as Earth rotates around the Sun [75]. 
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5.1.1. The dark halo model 


In order to compare the results of different experiments usually the Stan¬ 
dard Halo Model (SHM) is used. In the SHM the local DM density is 
Pdm — 0.3 GeV/cm 3 (its error is ±0.1GeV/cm 3 , see e.g. [76]). The lo¬ 
cal DM velocity distribution in the galactic reference frame /g(^g) is a 
Maxwell-Boltzmann distribution with zero average and dispersion Vq . trun¬ 
cated at the local escape speed v esc from our galaxy (see e.g. [77]). The 
WIMP velocity with respect to the Galaxy is vg = v + Vq + v E {t) where v 
is the WIMP velocity with respect to Earth, Vq is the velocity of the Sun 
with respect to the galaxy, and v E (t) is the velocity of Earth with respect 
to the Sun. Thus the velocity distribution with respect to the Earth is 

= f G (v + v® + v E (t)). (17) 


There are uncertainties in the parameters of the SHM (see e.g. [78]). Some 
usual values are Vq = 232 km/s [77], Vo = 220 km/s, and v esc = 544 
km/s [79]. The most recent Radial Velocity Experiment (RAVE) 2013 
results [80] give v esc = 533+| 4 km/s but values of v esc between 500 and 650 
km/s can be found in the literature. In the SHM the maximum average 
velocity of WIMPs with respect to Earth happens between the end of May 
and the beginning of June (see e.g. [78]). The minimum occurs six months 
later (although not exactly, because of the small ellipticity of Earth’s orbit). 

We expect the actual halo to deviate from this simplistic model. The 
local density and velocity distribution could actually be very different if 
Earth is within a DM clump, which is unlikely [81], or in a DM stream, or 
if there is a “Dark Disk [40, 39] (see section 2) in our galaxy. The DM of 
the Sagittarius Stream, tidally stripped from the Sagittarius Dwarf Galaxy, 
could be passing through the Solar system, perpendicularly to the galactic 
disk [83]. A large amount of DM clumps are expected to remain within 
the dark halo of our galaxy [81, 82], because haloes grow hierarchically, 
incorporating lumps and tidal streams from earlier phases of structure for¬ 
mation. However, clumps are more effectively destroyed by tidal effects 
near the center of the galaxy, thus most of them are far from the Sun. The 
chance that a random point close to the Sun is lying within a clump is 
smaller than 10 -4 [81]. “Debris flows”, which are spatially homogeneous 
structures in velocity, are expected from this complicated merger history, 
and they would also modify the velocity distribution primarily in the high 
velocity tail [84]. 
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5.1.2. Recoil energies and rates 

The maximum recoil energy of a target of mass Mt in an elastic collision 
with a WIMP of mass m is 


Emax = 2^V 2 /M T = q^ nax /2M T , (18) 


where \xt = toMt/(to + Mt) is the reduced mass. For light WIMPs with 
m <C Mt, Ht — m, and typically E max ~ 2keV (m/GeV) 2 (10 G eV/M t), 
since v ~ 10 -3 c. The threshold recoil energy in most detectors at present is 
O(keV), thus detectable WIMPs must have m > GeV. For heavy WIMPs 
with to Mt, M = Mt- The energy is large enough, E max = 2 At keV 
(we use Mt — At GeV if the nuclear mass number is At) but the limits 
die out because the WIMP flux decreases as 1/m. A DM candidate with 
sub-GeV mass, called “Light DM” (LDM), with m ~ MeV to GeV, could 
deposit enough energy through interactions with electrons (between 1 to 
10 eV) to be detected via electron ionization or excitation, or molecular 
dissociation [85, 86]. 

The typical momentum transfer in an elastic collision is q ~ /j,tV ~ 
O(MeV), which becomes q ~ MeV(m/GeV) for in <C Mt, and q ~ Aj-MeV 
for to ^ Mt- We note that 

q < l/I?Nucieu S ^ MeV (l60/Aj/ 3 ) , (19) 

thus WIMPs interact coherently with nuclei. I?Nucieus = 1-25 frn .4^/ 3 is the 
radius of a target nucleus (we recall that fm=10 _15 meters=(197 MeV) -1 ) 

-i/o 

and Aj, is a number of about 3 to 5 for most nuclei. If 1/q » I?Nucieus 
the nucleus interacts like a point-like particle. For larger q, the loss of 
complete coherence is taken into account by a nuclear form factor. For 
spin-independent (SI) interactions (see below) the form factor is the Fourier 
transform of the density of nucleons in the nucleus. A usual form for it is 
the Helm form factor for the distribution of charge (protons), assuming 
that the distributions of neutrons and protons are similar [87]. 

The recoil rate expected (in non-directional detectors) is given by the 
local WIMP flux nv times the number of target nuclei in the detector times 
the scattering cross section, integrated over the WIMP velocity distribution. 
In units of events/(unit mass of detector)/(keV of recoil energy)/day the 
expected differential recoil rate is 


dR v -v dRT 
dE R dE R 



daT 

dE R 


x nvf(v,t)d 3 v. 


( 20 ) 
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Here, Er is the recoil energy, T denotes each target nuclide (elements and 
isotopes), Ct is the mass fraction of nuclide T in the detector, Ct/M t is 
the number of targets T per unit mass of detector, da^/dER is the dif¬ 
ferential WIMP-nucleus scattering cross section and n = p/m. The local 
WIMP density p is the local DM density, p=Pdm , if the WIMP in question 
constitutes all the DM. If instead i?=fi W iMP /Cl DM <l, then p=Rp DM , as¬ 
suming that the DM close to Earth has the same composition as the DM 
as a whole. The local DM density and velocity distribution with respect to 
Earth, f(v,t), depend on the dark halo model adopted. 

The minimum speed u m i n a WIMP must have to communicate to the 
target T a recoil energy Er depends on the collision being elastic or in¬ 
elastic. In some particle models a DM particle of mass m may collide 
inelastically, producing a different state with mass m! = m + 5, while the 
elastic scattering is either forbidden or suppressed [88]). S can either be 
positive (“inelastic DM”, iDM, with “endothermic” scattering) [88] or neg¬ 
ative (“exothermic DM”, exoDM, with “exothermic” scattering) [89]. These 
types of interactions enhance the potential signal in some targets and sup¬ 
press it in some others: “endothermic” scattering favors heavier targets (so 
the iodine in DAM A/LIBRA is preferred over the germanium in CDMS) 
while “exothermic” scattering behaves in the opposite fashion (the silicon 
in CDMS is preferred over the xenon in LUX or XenonlOO). For elastic 
collisions 6 = 0. For <5| <C m, 



The differential recoil rate, Eq. (20), is not directly experimentally ac¬ 
cessible because of energy-dependent experimental efficiencies and energy 
resolution functions, and because what is often measured is a part E' of 
the recoil energy Er. The observable differential rate is 

§ = <e) f ll E x -£0 T (E K , E ')^. ( 22 ) 

Here E' is the detected energy, often quoted in keVee (keV electron- 
equivalent) or in photoelectrons. e(E') is a counting efficiency or cut accep¬ 
tance. Gr{Er, E') is a (target nuclide and detector dependent) effective 
energy resolution function that gives the probability that a recoil energy 
Er is measured as E'. It incorporates the mean value (E 1 ) = QrEr , which 
depends on an energy dependent “quenching factor” Qt(Er ), and the en¬ 
ergy resolution cte r (E'). These functions must be measured (although 
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sometimes the energy resolution is obtained from computation). 

WIMP interactions in crystals produce mostly phonons. Only a fraction 
Qt of the recoil energy goes on average into ionization or scintillation. For 
example, Qa e — 0-3, Qsi — 0.25, QiVa — 0.3, and Qi ~ 0.09. In noble gases 
such as Xe, a similar factor L e ff measures the scintillation efficiency of a 
WIMP relative to a photon. There are large experimental uncertainties in 
the determination of these parameters at low energies. 

We can see in Eqs. (22) and (20) that three elements enter into the 
observed rate in direct detection experiments: the detector response, the 
local dark halo model and, finally, the particle physics input, given through 
the cross section and mass of the DM candidate. 


5.1.3. The scattering cross section 

For contact interactions in the non-relativistic limit v —> 0 the differential 
cross section has the form 


dar (Tt(Er) , „ .. Mt 

= - = aT ^ R > 9 i ( 2 „ 2 - 

ClU/R t^max ZflrpV 


(23) 


for momentum transfer q and velocity-independent interaction operators. 
There are only two types of interactions of this kind, SI and SD. 

“Spin-independent” (SI) contact interactions are due to scalar or vec¬ 
tor couplings. The DM couples to the nuclear density. We can write 
a ? 1 (Eft)=crToF 2 (ER), where F 2 (Er) is the nuclear form factor, F' 2 (0)=1, 
(see e.g. [90]), and 

=\Z + (A- Z)(f n /f p )] 2 (^/fc 2 p )a p . (24) 

Here, a p is the WIMP-proton cross section while f n ,fp are the WIMP 
couplings to n,p. For (/ n // p )=l (isospin-conserving coupling), <jto = 
A 2 (Tt/t p) Up. Heavier nuclei are more neutron rich than lighter nuclei. 
Thus, the ratio {fn/fp) = —Zt/{At — Zt) cancels the coupling of the 
DM with the particular nuclide with atomic and mass numbers Zt and 
At, and changes the couplings to all other nuclei too (no choice can make 
the coupling with an element zero because of its isotopic composition). 
“Isospin-Violating” DM [91] with f n /f p = —0.7 minimizes the coupling 
with xenon (thus weakening the limits of XenonlO, XenonlOO and LUX, 
some of the strongest at present). f n /f p = —0.8 instead reduces maximally 
the coupling to germanium [92] (thus weakening preferentially the limits of 
CDMS and SuperCDMS). 
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Figure 4. Overview of 90%CL existing direct detection limits (solid lines) and future 
sensitivity regions for a WIMP with SI interactions and f p = f n constituting the whole 
of the DM, assuming the SHM, plotted in the (ra, cr p ) plane [74]. DAMA/LIBRA (light 
brown), CoGeNT (yellow), CRESST-II (pink) and CDMS-II-Si (light blue) signal regions 
for this WIMP, are shown. The level at which neutrinos of different sources would 
constitute a background to a DM signal, the “neutrino floor” (thick dashed orange lines) 
and regions of interest for particular DM candidates are also shown (see [74] for details). 


“Spin-dependent” (SD) interactions result from an axial-vector cou- 
pling. The DM couples to the nuclear spin density, leading to (see e.g. [90]) 

°t D (E r ) = 32 ^ 2 G 2 f [(J t + 1 )/Jt] [(S p )a p + (S n )a n ] 2 F$ d (Er). (25) 

Here, Jt is the nuclear spin, a p , n are the WIMP couplings to p and n and 
Esd(Er) is the nuclear form factor, with Fg D (0 )=1. (5 p , n ), the expecta¬ 
tion values of the p and n spin content in the target nucleus, are numbers < 
0 (1) that can differ easily by factors of 2 or more in different nuclear mod¬ 
els (see for example Ref. [93]). Since also the nuclear spins are of 0(1), SD 
cross sections are a factor Aj, smaller than SI cross sections. The bounds 
on the second are therefore typically better than the bounds on the first 
(Figs. 4, 5). 

There are many other types of possible DM-nucleus interactions be¬ 
sides the two mentioned, and many of them have been considered in recent 
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years to try to accommodate different hints of a DM signal in direct or 
indirect searches. All other interation operators contain extra powers of 
the momentum transfer q or WIMP velocity. Refs. [94, 95] list all possible 
operators in the non-relativistic limit up to 0(q 2 ). For example, a pseudo¬ 
scalar mediator yields a coupling (Sum • q) (Sn • q ), with two extra powers 
of q, dependent on the spin of the DM particle, Sdm, and the spin of the 
interacting nucleon, Sn- The mediators can be either heavy or light with 
respect to the momentum transfer. 

Starting from the fundamental interactions with which the DM particles 
couple to quarks, there are uncertainties on how to pass from quarks to 
protons and neutrons, and then to nuclei. Each type of interaction requires 
its own nuclear form factor, most of which are poorly known. 


5.1.4. The average inverse speed function r?(u m i„,t) 


With the contact differential cross-section of Eq.(23), Eq. (20) becomes 


dR v- c t t {e r ) r r f(v,t) 3 ' 

dE R ^ 2 m/4 P [J v>Vmin v 


E 


&t{E_ r) 

2 m/Xjn 


PV(V min, 


t), (26) 


where the last equality defines the velocity integral rj{v m \ n . t). Due to the 
revolution of the Earth around the Sun, rj(v m i n ,t) has an annual modu¬ 
lation generally well approximated by the first two terms of a harmonic 
series, 


Vaulin ,t) =’n 0 (v m i n )+1J 1 {v min )cOs(u}(t-to), (27) 

where to is the time of the maximum of the signal and oj = 27r/yr. In 
the SHM, r](v m \ n . t) is maximum when the average WIMP velocity with 
respect to Earth is maximum for z; m i„ > 200 km/s. For lower r; m i n values, 
rj(v m \n-1) is instead a minimum then. The time average rjo and the mod¬ 
ulated component of p enter respectively into the average and modulated 
parts of the rate in Eq. (26). Notice that the factor pr](v m i n ,t) includes all 
the dependence of the rate on the dark halo model for any detector. 


5.1.5. Hints and limits: Halo-Dependent and Independent analyses 

There are many direct DM detection experiments that are either running, 
in construction or in the stage of research and development (see e.g. [74] and 
references therein). They use different target materials and detection strate¬ 
gies (see Fig. 4). Three direct detection experiments, DAMA/LIBRA [96] 
(containing Nal), CoGeNT (Ge) [97] and CDMS-II-Si [98] have at present 
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Figure 5. Best bounds on WIMPs with spin-dependent (SD) interaction coupled only 
to protons, assuming the SHM. (Fig. from [74]). The IceCube limits come from the 
non-observation of neutrinos from the Sun (indirect detection), assuming annihilation 
into the mentioned channels. The others are direct detection limits. See [74] for details. 


claims of having observed potential signals of DM. DAMA/LIBRA finds in 
14 years of data an annual modulation at the 9.3<r C.L. with a phase com¬ 
patible with that expected from DM in the SHM. A few years ago CoGeNT 
found both an unexplained rate excess and an annual modulation (at the 2er 
initially, and later with a smaller C.L.), both attributable to WIMP inter¬ 
actions. CDMS-II (had Ge and Si) found 3 unexplained low-energy events 
in their Si component which could be due to DM. CRESST-II with an up¬ 
graded detector no longer finds an unexplained excess in their rate [99], as 
they had found in their previous 2010 results [100]. All other direct de¬ 
tection searches, including LUX (Xe), XENONIOO (Xe), XENONIO (Xe), 
CDMS-II-Ge, CDMSlite (Ge) and SuperCDMS (Ge), have produced only 
upper bounds on the interaction rate and annual modulation amplitude of 
a potential WIMP signal (see e.g [101] and references therein). 

It is thus essential to compare these data to decide if the potential DM 
signals are compatible with each other and with the upper bounds set by 
searches with negative results for any particular DM candidate. There 
are two ways of presenting and comparing direct detection data: “Halo- 
Dependent” and “Halo-Independent”. 

As we have seen, the rate observed in a particular direct-detection de¬ 
tector due to DM particles in the dark halo of our galaxy depends on 
three main elements: 1) the detector response to potential WIMP colli¬ 
sions within it; 2) the WIMP-nucleus cross section and WIMP mass; 3) the 
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local density p and velocity distribution /(v, t) of WIMPs passing through 
the detector. All three elements have considerable uncertainties. 

The “Halo-Dependent” data comparison method, used since the incep¬ 
tion of direct detection [102, 103], fixes the three aforementioned elements 
of the rate, usually assuming the SHM for the galactic halo, except for the 
WIMP mass m and a reference cross section parameter (7 re f extracted from 
the cross section (cr re f = cr p for SI interactions). Data are plotted in the 
(m,cr re f) plane (as in Figs. 4, 5 and the left panel of Fig. 6). All exper¬ 
imental collaborations use this method and assume WIMPs with SI with 
fn = fp to present their results. Fig. 4 shows that with these assumptions, 
the DAMA, CoGeNT and CDMS-II-Si regions almost overlap for “Light 
WIMPs” with m ~ 10 GeV, but they are all rejected by upper limits. One 
caveat is that in the “Halo-Dependent” plots it is usually assumed that the 
WIMP constitutes the whole of the DM. If it constitutes a fraction R of 
the DM, the vertical axis is Ra re f and not just cr re f- 

In the “Halo-Independent” data comparison method one fixes the el¬ 
ements 1) and 2) of the rate, except again for a reference cross section 
parameter cr re f, but does not make any assumption about the element 3), 
circumventing in this manner the uncertainties in our knowledge of the local 
characteristics of the dark halo of our galaxy [104, 105, 106] (see e.g. [107] 
and references therein). The main idea is that for a particular DM can¬ 
didate the interaction rate at one particular recoil energy depends for any 
experiment on one and the same function a le {pri(y m i n , t)/m = fj(v m i n ,t). 
Thus, all rate measurements and bounds can be mapped onto measure¬ 
ments of and bounds on the unique function fj(v m i n ,t) for a fixed WIMP 
mass m and plotted in the (v rrun , f]) plane. To be compatible, exper¬ 
iments must measure the same fj function. This method was initially 
developed for a SI WIMP-nucleus interaction [104, 105, 106] and only 
in [108] extended to any other type of WIMP-nucleus interaction. In 
this case the rate in an observed energy interval [E [, E 2 \ is written as 
= Je[ dE ' dR/dE’ = / 0 °° dv m i n TZ[ E ^ E ^(v min )fj(v in i n ,t) with a DM 
candidate and detector dependent response function TZ [106, 108]. 

The right panel of Fig. 6 shows the Halo-Independent comparison for 
the same candidate whose Halo-Dependent analysis assuming the SHM is 
shown in the left panel of the same figure: a WIMP with Isospin-Violating 
{fn/f P = —0.7) SI interactions and exothermic scattering (6 = —50 keV) 
(see [92] for details). The crosses in this panel represent potential rate and 
modulation amplitude measurements and we see that part of the three red 
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Figure 6. Halo-Dependent direct detection data comparison in the (m,<T p ) plane, in 
the SHM (left) and Halo-Independent data comparison in the (i> m i n , p<Tpt](v m in)/in) 
plane with m = 3.5 GeV (right) for an exothermic WIMP (<5 = —50 keV) with Isospin- 
Violating ( fn/fp = —0.7) SI coupling, for which the CDMS-II-Si signal (in pink-red) is 
compatible with all limits. The dashed gray lines in the right panel show the shape of 
the average (upper line) and the annual modulation amplitude of a p pri/m (lower line) 
in the SHM (with a v = 1 X 10 -40 cm 2 ). Fig. from [92] (see this reference for details). 


crosses corresponding to the three events observed by CDMS-II-Si escape 
all bounds. This is consistent with the analysis in the left panel, where part 
of the CDMS-II-Si region is also seen to escape all bounds. 

So far the results have indicated that one of the potential direct DM 
detection signals at a time could be compatible with all negative results 
for particular DM candidates [95, 109, 92], but not all of them. For 
example, a signal in DAM A/LIBRA is favored by a magnetic dipole- 
moment coupling [95] or a spin-dependent coupling to protons, such as 
(Sdm ■ <l){Sp ’ q) [109]. Both favor couplings to Na and I and disfavor 
couplings with Xe and Ge. “Magnetic Inelastic DM” [95] may work as a 
candidate for DAM A/LIBRA because inelasticity further disfavors lighter 
targets (better I than Ge). 

The situation will be clarified with more data and very fast advances are 
expected in direct detection in the next decade. Fig. 4 presents an overview 
of existing direct detection limits (solid lines) and future sensitivity regions. 
The next generation of multi-ton experiments will reach the rate level at 
which neutrinos of different sources (atmospheric, the Sun, the diffuse back¬ 
ground from supernovas or DSNB) will constitute a background to a DM 
signal, called the “neutrino floor” (see [74], for details and references). 
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5.2. Indirect Detection of WIMPs 


The DM annihilation or decay products searched for are high energy neu¬ 
trinos produced by WIMPs captured inside the Sun or the Earth, photons 
and anomalous cosmic rays, such as positrons and antiprotons, which do 
not come from astrophysical sources. Indirect searches assume that DM 
particles annihilate or decay. Here we assume that they annihilate (e.g. 
because they are self-conjugated Majorana fermions or bosons). 


5.2.1. Indirect detection through neutrinos from the Sun 

As the Sun and Earth move through the dark halo of our galaxy, WIMPs 
can occasionally scatter with the baryonic matter in one of these bodies and 
loose enough energy to become gravitationally bound to it. The capture 
rate Tc is proportional to the scattering cross section and the local DM 
number density, Tc ~ asn. As the number of captured WIMPs N = (T c t) 
increases with time, the annihilation rate increases too as Ta ~ oaN 2 , 
and the number of WIMPs inside the Sun or Earth changes as dN/dt = 
Tc — 2Ta- The capture rate is constant, since the captured number of DM 
particles is too small to affect the DM density and distribution in the halo. 
Given enough time, the annihilation rate grows enough to compensate the 
capture, Ta = Tc/2, so dN/dt = 0 and a constant equilibrium number is 
reached. For typical WIMPs, equilibrium is reached in the Sun within the 
lifetime of the Sun (4 x 10 9 y) but not in the Earth [110, 90]. 

From the annihilation products, only neutrinos can escape from the Sun 
or Earth and travel to large-volume underwater or ice “neutrino telescopes”, 
where they could be detected. These are neutrinos with energies much 
larger than those otherwise produced in these bodies. Besides, a signal from 
the Sun would follow its position, thus eliminating all possible backgrounds. 

The best limits on neutrinos from the Sun come from the IceCube tele¬ 
scope [111], an instrumented kilometer cube of ice in Antarctica (with a 
more densely instrumented core, called DeepCore). A smaller underwater 
neutrino telescope in the Mediterranean, Antares [112], is a prototype for a 
possible “KiloMeter-cube Neutrino Telescope”, KM3NeT, to be built some¬ 
where in the Mediterranean. In the limit of equilibrated capture and anni¬ 
hilation, the flux of neutrinos only depends on the rate at which WIMPs are 
captured. This rate depends on the scattering cross section. Since direct 
detection depends on the same quantity both types of bounds can be com¬ 
pared, as in Fig. 5. Because the Sun contains primarily hydrogen, the limits 
on the spin-independent (SI) scattering cross section, Eq. (24), obtained by 
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IceCube are not competitive with those obtained in direct detection experi¬ 
ments. Spin-dependent (SD) WIMP-nucleus cross sections, Eq. (25), are of 
the same order of magnitude for hydrogen and heavy nuclei, and therefore 
the Ice Cube limits on these cross sections are very important, as shown in 
Fig. 5. 

The indirect limits due to neutrinos in the Sun are particularly impor¬ 
tant for Self-Interacting DM (defined in section 3) because the capture rate 
is enhanced by the scattering of halo DM with captured DM, once the 
number of captured particles becomes large enough [113]. 


5.2.2. Indirect detection through photons 

Gamma-ray astronomy is done with ground and space instruments. The 
Fermi Space Telescope was launched in 2008. Its main instrument is the 
Large Area Telescope (LAT) which detects photons between 20 MeV and 
300 GeV. Photons with energy above 20 GeV and up to several TeV 
are detected by ground-based Air Cherenkov Telescopes (ACT): HESS in 
Namibia, MAGIC in Las Palmas and Veritas in the US. On the planning 
stage is a large array of ACTs, the CTA (Large Telescope Array), which 
could detect photons from 10’s of GeV to above 100 TeV. 

Photons reveal the spatial distribution of their sources. The Universe is 
totally transparent to photons below 100 GeV. At higher energies photons 
interact with infrared and optical backgrounds, but below 10’s of TeV still 
arrive to us from cosmological distances (at higher energies they are ab¬ 
sorbed by the CMB; only above 10 10 GeV, the range of Ultra-High Energy 
Cosmic Rays, can they again reach us from cosmological distances). 

Monochromatic photons can be produced through \X ~► 77 (or 7 Z) 
with energy equal (or close) to the annihilating WIMP x mass m. Detec¬ 
tion of this monoenergetic spectral line would be a “smoking gun” signa¬ 
ture. However, usually (but not always) these processes happen only at the 
one-loop order and they are suppressed with branching ratios 10 ~ 3 — 1CU 4 . 
Secondary photons, in particular from pion 7 r° —► 77 decay, would be pro¬ 
duced with a spectrum whose cutoff at high energies is m. So a cutoff in 
the observed spectrum would be a signature of DM annihilation, but not 
exclusively. The spectrum of astrophysical sources also has a cutoff, be¬ 
cause there is a limit to the energy with which emitted particles can be 
produced. 

The 7 -ray flux (number per unit area, time and energy) expected from 
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Figure 7. Radial DM density distribution for different galactic halo models. The central 
density value is fixed by the requirement of the rotation speed of the Sun [115] 


the annihilation of DM particles coming from a particular direction is 

<J> 7 (£ 7 ) ~ (a A v) ^ f ^dl(8)dd , (28) 

U-tyj ./line of sight 

where (a A v) is the annihilation cross section times relative speed at the 
source, dN 1 /dE 1 is the 7 -ray spectrum per annihilation (for example, it 
would be a delta function for XX ► 77 ). The integration of the DM number 
density squared as function of distance, (p(r)/m) 2 , is along the line of sight 
and over the angular aperture or resolution of the detector. 

Since the annihilation rate depends on the square of the DM density, 
high density regions such as the galactic center (GC), DM clumps, dwarf 
galaxies and other galaxies and galaxy clusters would boost the rate. Thus 
a signal is expected from them. 

The amount of DM in the GC is uncertain (see Fig. 7). In fits to nu¬ 
merical simulations of structure formation containing CDM only, Navarro, 
Frenk and White (NWF) found a universal density profile [114], at all scales 
(large clusters, galaxies and dwarf galaxies) which goes as 1 /r 3 at large dis¬ 
tances from the GC and as 1/r close it. This is called a “cuspy profile.” 
It is expected to be modified close to the center by the effect of baryons, 
which are dissipative, and thus fall to the center and dominate over the DM. 
Fits to observed rotation curves of galaxies find instead “cored profiles”, 
in which the DM density becomes constant close to the GC. A commonly 
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used halo profile is the “generalized NFR”, 

p(r) = Po {r/rsY 1 [1 + (r/r s )] 7 ~ 3 , (29) 

where r s ~ 20 kpc is called the “scale radius” and the central density value 
Po is fixed by imposing that the DM density at the distance of the Sun from 
the GC is equal to the value inferred by observations (~ 0.3 GeV/cm 3 ). For 
r <C r Sl p ~ 1/r 7 , while for r r s , p ~ 1/r 3 . The NFW profile has 7 = 1 
and a cored profile has 7 ~ 0. When considering the region of the GC, the 
most important feature of the halo profile is the inner slope, 7 . As shown 
in Fig. 7 [115], in some halo profiles, this slope can be considerably steeper 
than in the NFW profile. E.g. the Moore profile has 7 = 1.5 [116]. 

Several hints for DM have been found at the GC. We review these 
potential signals below. 

5.2.3. The INTEGRAL 511 keV signal 

The satellite INTEGRAL, launched in 2002, confirmed the emission of 511 
keV photons from the GC, a 40 year old signal first observed by balloon- 
born 7 -ray detectors. It is clearly due to a non-relativistic e + annihilating 
with an e~ almost at rest. Initially the emitting region seemed spherically 
symmetric, as expected for a DM signal, but in 2008 INTEGRAL [117] 
revealed that is it not, and found evidence of a population of binary stars 
consistent with being the e + source. These observations decreased the mo¬ 
tivation to consider DM as the origin of the signal. Special DM candidates 
were proposed to explain this signal: “Light DM”, LDM, with MeV mass, 
annihilating as XX e + e _ into e + e _ almost at rest [118]; “eXciting DM”, 
XDM, in which a DM particle x! decays into a lighter one x with mass dif¬ 
ference ml — m = 5 ~ MeV and an e + e _ almost at rest, x ’ —► X e+e ~ [119]. 

5.2.4. The “WMAP Haze” 

This is an excess of microwave emission in the inner 20 degrees, about 1 
kpc, around the GC. It was discovered by D. Finkbeiner in 2003 [120]. It 
can be explained as synchrotron radiation from e~ and e + (accelerated in 
magnetic fields) produced in astrophysical sources or by the annihilation of 
DM particles. Almost any WIMP annihilating into e~ or e + could produce 
it. However, it is now considered part of the “Fermi Bubbles” (discovered 
in 2009 also by Finkbeiner and collaborators) [121]. These are two large 
structures of 7 -ray emission of 8 kpc in diameter each, extending to both 


33 




Figure 8 . Spectrum of the extended GeV excess from the GC (red crosses) and several 
fits with DM annihilating into bb and t+t~. Fig. from [124]. Notice the required inner 
slope of the galactic density profile 7 = 1.2 (a cuspy profile) necessary to fit the data. 


sides of the galactic plane. They could be due to an early period of strong 
jet-like activity of the black hole at the GC, which is now dormant. 


5.2.5. The extended GeV excess 

In 2011, Hooper and Linden [122] found by subtracting from the Fermi data 
all known contributions, an unexplained extended excess of GeV photons, 
peaking at 2 -3 GeV, coming from the GC (see Fig. 8 ). The existence of 
the excess was confirmed by several other groups [123, 124] and later also 
found in the inner galaxy (up to 10 degrees from the GC) [125]. 

This signal, shown in Fig. 8 , can be interpreted as possible evidence of 
DM particles with mass of 7-12 GeV annihilating to t + t~ (among other 
leptons) or with mass of 22-45 GeV annihilating to quarks, with an annihi¬ 
lation cross section close to the value ~ 10 -26 cm 3 /s required by thermal 
WIMPs at decoupling (see Fig. 9). In both cases to fit the observed signal 
the galactic profile must be “cuspy” with an inner slope 7 ~ 1.2 to 1.4. A 
cuspy DM halo profile enhances the rate from the GC considerably with 
respect to a core profile. Fig. 7 shows that the ratio of the central densities 
of the “cuspy” Moore (7 = 1.5) to the “cored” 7 = 0 profiles is about 30, 
which translates into a factor of 900 in the rate predicted from the GC. 

The GeV excess could also be explained by unresolved millisecond pul- 


34 








3 


-bb, x 2 = 26.4 

-cc, X 2 =23.7 



m x (GeV) 


Figure 9. DM mass and cross section required to explain the extended GeV emis¬ 
sion from the GC if the DM annihilates into the shown quark-antiquark modes. Fig. 
from [125]. Notice the required inner slope of the galactic density profile, 7 = 1.26. 


sars [127, 123, 124, 126]. These are rotating neutron stars which have been 
spun-up, through accretion from a companion star with which it forms a 
close binary system, to a period of 1-10 milliseconds. They emit X-rays 
and 7 -rays, possibly from the matter being accreted. Their distribution in 
or near the GC is not known, as they cannot be observed. There are also 
uncertainties in the photon spectrum they emit. Several assumptions are 
needed to fit the observed GeV excess. 

5.2.6. The 130 GeV line 

A line at this energy with 3.5cr significance was found in 2012 in the Fer- 
miLAT data by C. Weniger and collaborators [128] coming from the GC. 
This monochromatic line could have been produced by DM annihilation 
into 77 . A hint of a second line was found at 111 GeV [129] which could 
have corresponded to annihilation of the same particles into Z° 7 . However, 
no evidence of the line was found elsewhere where it was expected if due 
to DM annihilation, e.g. in dwarf galaxies [130]. The line was also found 
where it could not be if due to DM annihilation: in Earth’s Limb [131] and 
in the vicinity of the Sun [132] (Earth’s Limb is the edge of the disk of the 
Earth, where 7 -rays can be produced by cosmic rays). The last two are 
clear indications of a detector effect. The FermiLAT paper [133] about the 
line finds a hint of it at the 1 . 6 er level from the GC, and also at 2<r in the 
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Figure 10. Map of known dwarf galaxies overlaid on Fermi data. Fig. from [136]. 

Limb. Another troubling indication for a DM interpretation is that the ac¬ 
cumulated significance of a real annihilation line should increase with time 
but it is decreasing for this signal [134]. Thus, the line seems spurious, and 
possibly due to an experimental effect, although the issue is not entirely 
resolved yet [135]. 

5.2.7. FermiLAT limit on thermal WIMPs from dwarf galaxies 

The GC is a complicated place, with large uncertainties in the DM profile, 
and other powerful photon sources. Other overdense regions might provide 
a cleaner signal. A large amount of DM clumps are expected to remain 
within the dark halo of our galaxy. CDM simulations of structure formation 
(the Via Lactea II [82] and the Aquarius Project [81]) suggest that 0(10) 
clumps could be discovered by Fermi, although none has been observed so 
far. Dwarf galaxies are simpler sources because they are the most DM- 
dominated structures observed so far. A map of the known dwarf galaxies 
overlaid with the Fermi gamma-ray map is shown in Fig. 10. 

The best limit on thermal WIMPs from the Fermi-LAT collaboration 
comes from 15 stacked dwarf galaxy images [136]. When observations on 
each of several very similar objects are not good enough, stacking their 
images increases the statistics and allows one to obtain better limits. Fig. 11 
shows the upper limits on (gav) assuming exclusive annihilation into the 
pairs of leptons or quarks indicated in each panel. To have an acceptable 
density < Cl DM , thermal WIMPs must have (gav) > 3 x 10 -26 cm 3 /s at 
decoupling. This requires the limit to be above the horizontal dashed line 
in Fig. 11. Both limits together impose m > 20 GeV for thermal WIMPs 
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(from the t + t~ and uu modes) and a limit m > 0(100) GeV is expected 
in 4 y [136]. 

The caveat to this limit is that (< 7aV ) at decoupling and now may not be 
the same, because the characteristic WIMP speed at decoupling is v ~ c/3 
and in the galaxy is v ~ 10~ 3 c. For s-wave annihilation (av) is indepen¬ 
dent of v and the limit applies, but for p-wave annihilation (av) ~ v , and 
the limit does not apply. Also, there could be a “Sommerfelcl enhance¬ 
ment” of the present annihilation which would invalidate the limit. This 
enhancement is due to a modification of the wave function of annihilat¬ 
ing low velocity particles due to attractive long-range forces. The particles 
form an almost bound state which enhances the annihilation cross section. 


5.2.8. Anomalous cosmic rays 

Positrons and antiprotons would be produced in DM annihilations in equal 
numbers as electrons and protons. They are an interesting potential signal 
of WIMP annihilation because there is not much antimatter in the Universe. 
Unlike photons, which travel in straight lines and are not absorbed for 
energies < 100 TeV, e + and e~ interact with the magnetic fields of the 
galaxy. They rapidly (within a few kpc) loose energy through the emission 
of synchrotron radiation and inverse Compton scattering interactions with 
photons (in Compton scattering a high energy 7 interacts with an e almost 
at rest; in inverse Compton a high energy e interacts with a 7 with much 
smaller energy and produces a lower energy e and a higher energy 7 ). p 
and p suffer convective mixing and spallation. They propagate further than 
electrons, but still only from a fraction of the size of the galaxy. 

Balloon-born experiments detecting positrons have found since the 
1980’s a possible excess over secondary cosmic ray fluxes. It was called the 
“HEAT excess” and is now called the “PAMELA excess”. PAMELA (Pay- 
load for Antimatter Matter Exploration and Light-nuclei Astrophysics), a 
satellite carrying a magnetic spectrometer that was in operation from 2006 
to 2011 , reported in 2008 an excess in the positron fraction e + /(e + + e~) 
in the 10 to 100 GeV energy range [137] compatible with the HEAT ex¬ 
cess. The result was confirmed by FermiLAT [138] and more recently by 
AMS-02 [139] (Alpha Magnetic Spectrometer), a cosmic ray research mod¬ 
ule mounted to the exterior of the International Space Station, in operation 
since 2011. Fermi does not have a magnet to distinguish positively and neg¬ 
atively charged particles, but they cleverly used the magnetic field of the 
Earth [138], 
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Figure 11. Fermi upper limits from 15 stacked dwarf galaxy images on the DM on 
{<jav), assuming exclusive annihilation in each of the modes shown [136]. The lower 
limit on (<tav) at freeze-out for thermal WIMPs is also shown (horizontal dashed line). 


The rapid rise shown in Fig. 12 in the cosmic ray e + fraction above 10 
GeV measured by PAMELA and AMS indicates the existence of nearby 
primary sources of high energy positrons, such as pulsars or annihilat¬ 
ing/decaying DM. The spectrum of secondary positrons produced through 
the collisions of cosmic rays in the interstellar medium is predicted to fall 
rapidly with energy, and thus is unable to account for the observed rise. 

It had been proposed that secondary positrons could be produced and 
then accelerated in nearby supernova remnants, potentially explaining the 
observed rise without the need of primary sources. If so, other secondary 
cosmic ray species (such as boron nuclei and p) should also be accelerated, 
leading to rises in the boron-to-carbon and p/p ratios. Recent measure¬ 
ments show no sign of such a rise, what disproves this mechanism [140]. 

It is very important for the DM interpretation that no excess in p has 
been found [141]. This means that the DM particles should be “leptophilic”, 
i.e. only annihilate (or decay) into leptons and not quarks, either because 
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Figure 12. The positron fraction e+/(e + + e ) measured by AMS (red), PAMELA 
(blue) and Fermi-LAT (green). Fig. from [139]. 


the DM carries lepton number or because of kinematics. 

No endpoint of the e + excess was observed by PAMELA or AMS, which 
would be an indication of the mass m of the DM particle. HESS measured 
the electron plus positron (which they call generically electrons) spectrum 
showing a soft decrease with energy from 1 to 10 TeV [142]. Together these 
observations require DM with m — few TeV annihilating not directly into 
e + e _ pairs (which would generate a sharp spectrum) but into leptons of 
the 2nd or 3rd generation: t + t~, 4/z, 4r or into pions [143]. This could 
be explained if the DM annihilates into a pair of bosons (j> with < 1 
GeV, XX < t ,( t > i and then <f> decays primarily into or pions because 

of kinematic reasons [37]. Due to the high multiplicity of such processes, 
the resulting e + energy distribution at injection is soft, as needed. 

Since (j> is so light, it must be very weakly coupled to ordinary matter 
to have escaped detection. The paper of Arkani-Hamed et al. titled “A 
Theory of Dark Matter” [37] containing this idea proposed the existence of 
a complicated dark sector in which the particles <fi are the light gauge bosons 
(“dark photons”) of a hidden gauge symmetry. In this model WIMPs with 
mass m ~500 to 800 GeV have exited states with small mass differences <5 
between 0.1 to 1 MeV, which can be “inelastic DM” (see section 5.1 and 
Eq. (21)) and “eXiting DM” (see the INTEGRAL signal above). 

As if all these necessary properties the DM should have to account for 
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the PAMELA excess are not complicated enough, the annihilation rate 
must be larger than expected for thermal WIMPs by a boost factor B ~ 
10 to 10 3 . The boost factor B is simply a factor that must multiply the 
spectrum obtained from the annihilation of thermal WIMPs and that must 
somehow be explained. Astroplrysical enhancements due to nearby over- 
dense structures of DM cannot be larger than a few. Most DM clumps are 
expected to be far from the Sun [81]. A possibility is that the annihila¬ 
tion cross section is larger in the dark halo at present, but that it had the 
value necessary for thermal WIMPs at freeze-out. We already mentioned 
this idea of a Sommerfeld enhancement of the present annihilation. The 
exchange of a light boson <f) may produce a Yukawa potential and as the 
annihilating particles approach at very low relative velocity they almost 
form a bound state, which enhances the cross section [37]. No boost factor 
at all is needed if WIMPs have a large annihilation cross section in both the 
early Universe and the dark galactic halo near Earth, which would produce 
a too-small relic abundance for thermal WIMPs but could be fine if the 
pre-BBN cosmology is non-standard [144]. 

Besides all these requirements on the type of DM, constraints imposed 
by its annihilation in the GC are only compatible with halo models that 
predict a relatively small amount of DM in the GC (cored profiles) [145]. 
Decaying DM has also been considered (see e.g. [145, 143]). It must decay 
mostly into leptons of the 2nd or 3rd generation, have multi-TeV mass and 
a very long lifetime, r ~ 10 26 s. 

Upcoming AMS data may help to settle the origin of the positron excess 
not only by increasing statistics and extending studies to higher energies, 
but also by further constraining any anisotropy in the positron and electron 
flux. If the origin of the positrons is one of the pulsars nearby, there should 
be an anisotropy at some level. 


5 . 3 . WIMP searches at the LHC 

DM particles escape detection at colliders, thus they are characterized by 
missing transverse energy in collider events. One way to search for DM 
particles at colliders is to know the particle spectrum and how the DM 
couples to other charged and/or strongly interacting particles that can be 
directly observed in the detectors. This is how complicated decay chains 
predicted in supersymmetric models are searched for. Another way is to 
search for the production of a pair XX of WIMPs and one visible particle 
emitted either by the initial or the intermediate SM particles to detect the 
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Figure 13. Monojet diagram (from [151]). Table 1. Some Effective Field Theory op¬ 
erators [152] coupling a WIMP x to quarks q or gluons g , valid when the energy of the 
process is smaller than M*. If the process is due to a mediator of mass M and coupling 
A, the operator represents a contact interaction with propagator A 2 /M 2 = 1/Mj. 


event [146, 147, 148]. If the one observable particle emitted is a photon, 
it is a “monophoton” event; if it is a gluon, it is a “monojet” event (see 
the diagram in Fig. 13). Other mono-particle events such as mono-W’s 
(leptons), mono-Z’s (dileptons), or even mono-Higgses, have been consid¬ 
ered [149, 150], 

So far the limits obtained using monophoton, monojet, or other mono¬ 
particle events have been derived assuming effective couplings. Non- 
renormalizable operators from effective field-theory analysis such as those 
in the Table 1 [152] have been used to couple the DM particle \ to quarks or 
gluons. These effective couplings are valid when the energy of the process 
is small compared with the scale M*. If a process is due to the exchange of 
a mediator of mass M and coupling A, the mediator is assumed to be heavy 
with respect to the momentum transfer q of the LHC partonic reaction, 
M » q , so that its propagator reduces to a constant A 2 /M 2 = 1/M*. 
This approach is limited because it does not include possible interference 
between different operators, and it does not deal with lighter mediators. 

The same effective couplings in Table 1 can be used to compute the 
WIMP interaction with nuclei in direct searches, and the results can be 
compared in plots like the one shown in Fig. 14 (from [148]; see this reference 
for details). This type of plot must be understood with care. If the WIMP 
in question constitutes a fraction R of the local DM density, the direct 
detection limits must be multiplied by 1 /R, while the LHC limits do not 
change. More importantly, while it is valid to include the direct detection 
limits when presenting the LHC limits derived from contact interactions, 
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the reverse is not correct. The reason is that mediators that are heavy with 
respect to typical LHC partonic energies, M >100’s GeV, are also heavy in 
comparison with the typical MeV momentum transfer in direct detection 
experiments. The opposite is not true: what is a contact interaction in a 
direct detection experiment may very well not be such at the LHC. If the 
mediator is light enough, it could itself be produced at the LHC or other 
colliders, and the analysis of the collider data would be different. 

Once the mediator is light, its couplings with different SM particles 
must be specified to study it. The problem is then that many theories 
will show a common low energy behavior when the mediating particles are 
heavy compared to the energies involved. Each effective contact interaction 
corresponds to many different possible particle models for the mediator. 
Without going to a complete theory, like a complete supersymmetric model, 
an intermediate step is to try to study simplified models for DM at the LHC. 
These incorporate all the known constraints on different interactions, and 
classify mediators according to whether they propagate in the s-channel 
or the t-channel, or by the way the DM relic density occurs [153, 154]. 
Significant work remains to be done in this direction. 

5 . 4 . Complementarity of WIMP Searches 

Direct, indirect and collider DM searches are independent and comple¬ 
mentary. They differ in essential characteristics and rely on different DM 
properties to see a signal. If a compelling DM signal is discovered, comple¬ 
mentary experiments will be necessary to verify the initial discovery, and 
to determine the actual abundance in the Universe and properties of the 
particle in question. 

Attention must be paid to the assumptions that go into the quantitative 
comparison of limits coining from the three types of searches, which is 
necessarily model dependent. For example, as part of the Snowmass 2013 
study a simple contact interaction model was assumed [155]: 

Jp X7 m 75X J2 «7m759 + ^ XXG a ^G“„ + Y xi^X Y hn* > (30) 

where the interactions with quarks mediate spin-dependent direct signals, 
and those with gluons mediate spin-independent direct signals. The coeffi¬ 
cients M q . M g , and M\ characterize the strength of the interaction with the 
respective SM particle, and were chosen so that the combined annihilation 
cross section into all three channels provides the correct relic density for 
the thermal WIMP \ to constitute the whole of the DM. With this simple 
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Figure 14. LHC bounds (monojet limits and others) on SI DM-nucleon scattering com¬ 
pared to direct detection limits [148]. For DAMA and CoGeNT 90% and 3cr contours 
and for CRESST la and 2 a contours are also shown. See [148] for details. 


DM interacting with quarks 



Figure 15. DM discovery prospects for current and future direct detection, indirect 
detection and particle colliders on the cross section a(xx quarks)/ ct^ as function of 
the WIMP mass m x for WIMPs with the simple contact interaction in Eq. (30) (see [155] 
for details). Here a t ^ = 3 x 10 -26 cm 3 /s taken as the reference value for the annihilation 
cross section, is the value required for a thermal WIMP to account for all of the DM. 

model it is possible to compare the limits imposed by direct, indirect and 
collider searches as well as their reach (see Fig. 15 from [155]). 
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Figure 16. Results from a scan of the 19 pMSSM parameters plotted in the mass m x vs. 
R&si plane. R=tt x /Q,DM is the density fraction of neutralinos x in the DM, and a si is 
the WIMP-proton cross section for SI WIMP nucleus interactions. Each point is a model: 
those in the reach of direct detection (the future LZ, whose discovery limit is shown with 
a black line) are green, those within the reach of indirect detection (FermiLAT and the 
future CTA) are blue, those in the reach of both are red. The magenta points are models 
tested by LHC searches which escape detection in direct or indirect detection. The gray 
points escape all searches in the near future. See [156] for details. 



— LZ SI exclusion, 1000 Days • Excluded by ID but not DD 

• Survives DD, ID, and LHC • Excluded by DD and ID 

• Survives DD and ID but not LHC • Excluded by DD but not ID 


mix") (GeV) 


The limits in Fig. 15 still depend on the particular halo model and the 
details of the detector response assumed (see [155] and references therein 
for details). Besides, the limits would become very different if the particular 
DM candidate studied would constitute only a fraction R of the DM. In 
this case, the upper limits from direct detection on the scattering cross 
section as become larger by a factor of 1/R (since the rate is ~7?<Tg), the 
indirect detection limits on the annihilation a a would change by this factor 
squared, 1/7? 2 , and the LHC limits would remain unchanged. 

The complementarity of searches has also been studied in a complete 
model, the pMSSM (phenomenological Minimal Supersymmetric Standard 
Model), scanning over 19 parameters for DM mass values of 50 GeV to 4 
TeV (see Fig. 16 from [156]). The DM candidate is the lightest neutralino 
Xi, which is its own anti-particle. 
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In Fig. 16 each point represents one particular supersymmetric model; 
there are 200,000 points. Within each model, the DM interactions are 
completely specified and all relevant DM signals can be computed. Models 
which produce a subdominant component of neutralino DM have also been 
included. The vertical axis of Fig. 16 shows the product of the computed 
density fraction of neutralino DM, R=Q x /£Idm, times the cross section. 
The models are divided into categories indicated by the color of each point, 
depending on whether the DM candidate is within the reach of future di¬ 
rect detection experiments (green points) such as LZ (the future upgrade 
of the LUX experiment, whose reach is indicated with a black line), or 
within the reach of indirect detection experiments (blue points), in partic¬ 
ular FermiLAT and the future CTA, or both (red points). The magenta 
points represent models that escape detection in the direct and indirect 
DM experiments just mentioned, and are tested at the LHC. Fig. 16 shows 
clearly that the three different DM probes are necessary to test most of 
the supersymmetry models in this scan. It is interesting to see that the 
DM candidates in the gray models escape all searches in the near future 
(see [156] for details.) 


6. Axions as DM candidates 

In the SM there is a source of CP violation besides the phase of the 
Kobayashi-Maskawa mixing matrix. It is associated with the so called 8 pa¬ 
rameter, which provides the magnitude of a CP-violating interaction among 
gluons allowed by all general principles. Experimental bounds on electric 
dipole moments constrain this parameter to be very small, 8 < 1CD 10 , while 
a priori we would expect 8 ~ 1. This unexplained smallness of the 8 pa¬ 
rameter constitutes the so-called “strong CP problem”. The only known 
viable solution of this problem is the Peccei-Quinn mechanism. It consists 
of augmenting the SM (e.g. by the addition of a second doublet Higgs 
field, but there are other ways) so that the resulting Lagrangian has a new 
global U(1)pq chiral symmetry [157]. This global symmetry must be spon¬ 
taneously broken, and this generates a Goldstone boson, the “axion” [158]. 
The Goldstone mode is the component a of the scalar field along the degen¬ 
erate orbit of minima of the characteristic inverted mexican hat potential. 

In order for axions to be phenomenologically viable the spontaneous 
symmetry breaking (SSB) scale, f a , must be large, f a 100 GeV. QCD 
effects violate explicitly the PQ symmetry and generate a potential m a a 2 / 2 
for the axion field a = 8f a with m a ~ ^qcd/ fai whose minimum is at 
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Figure 17. Excluded ranges and regions of interest in terms of axion mass m a or PQ 
spontaneous symmetry breaking scale f a oc l/m a . The mass range //eV to meV cor¬ 
respond to the so called “classic region” of CDM axions and the mass range below 
(larger values of f a ) is the “anthropic region” of CDM axions. CAST and ADMX are 
experiments that search for axions. Fig taken from [162], See also [160]. 


9 = 0. The small explicit breaking produces a tilt in the orbit of the minima, 
which has then a minimum value and a non-zero curvature close to it. The 
PQ symmetry therefore solves the strong CP problem by transforming the 
9 parameter into a field that has a minimum at 9= 0, simultaneously giving 
a small mass to the axion, which becomes a pseudo-Goldstone boson. For 
good reviews on axions with complete reference lists see [159, 160]. 

An initial displacement cq = f a 9.i of the axion held from its minimum 
will result in coherent oscillations of this held with frequency m a c 2 /h. The 
present energy density in these oscillations is [161] (recall m a ~ 1 / f a ) 

1.184 /1n 1.184 


n n h 2 = 0.195 Of 


fa 


10 12 GeV 


= 0.105 9 2 


10 /neV 


(31) 


and behaves as CDM. Other types of light bosons, called “Axion-Like Par¬ 
ticles” (ALPs), are pseudo-Goldstone bosons of other broken global sym¬ 
metries which do not couple to QCD, and can also be good DM candidates. 
They can acquire masses through their interactions with strongly-coupled 
hidden sectors or explicit breaking of the associated global symmetry [162]. 
ALPs together with very light “hidden” gauge bosons (“dark photons”) are 
generically called “Weakly Interacting Slim Particles” (WISPs) [163]. 

Axions can also be HDM for “large” masses m a ~ eV (when produced 
thermally via a coupling with pions such as amm). The CERN Axion 
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Figure 18. CDM axion parameter space, using the standard computations of the axion 
relic abundance [164]. The yellow region is excluded by astrophysical and CMB con¬ 
straints. The now rejected BICEP2 claim would have implied a Hubble scale during 
inflation indicated by the green vertical band. Fig. from [164]. 


Solar Telescope (CAST) experiment looks for this type of axion, which 
could be emitted by the Sun. Axions as CDM must have much smaller 
masses. They are searched for with resonant cavities in the Axion DM 
experiment (ADMX) through the (model dependent) axion coupling with 
photons. Fig. 17 [162]. shows the excluded ranges and regions of interest in 
axion mass m a and PQ SSB scale f a (in particular the ranges within reach 
of CAST and ADMX). 

The mechanism of production of CDM axions depends on when the 
PQ SSB occurs relative to inflation: after or before. When the PQ SSB 
happens, the position along the orbit of minima where the field ends up 
is different in causally disconnected volumes in the Universe. This is the 
initial value cq = f a 0i of the field oscillations which happen much later, 
after the explicit PQ symmetry breaking becomes important. 

If inflation happens after the PQ SSB, our whole Universe comes from 
only one causally connected volume before inflation. The so-called “mis¬ 
alignment” Oi is then the same in the whole visible Universe. For this to 
happen the condition is f a > Hj/2n, where the Hubble expansion Hj is 
the only dimensional relevant quantity during inflation (see e.g. [164]). In 
this case, if Oi — 0(1), axions would have the DM density for m a ~ 10/reV 
and f a ~ 10 12 GeV (see Eq. (31)). This is called the “classic region” (see 
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Fig. 17). If instead 0* <C 1, Eq. (31) shows that axions should have a much 
smaller mass (and larger f a ) to be the DM. This is called the “anthropic 
region” (see Fig. 17). 

With the standard axion production mechanism, the possibility of f a > 
Hj/2tt would have been excluded [164] by the recent BICEP2 [12] claim 
(now rejected [11]) of a measurement of gravity waves produced during 
inflation, as shown in Fig. 18. The BICEP2 claimed measurement would 
have fixed the Hubble expansion during inflation to be Hj ~ 6 x 10 14 
GeV, for which limits on axion isocurvature fluctuations in the Universe 
forbid f a > Hj/2n. Isocurvature fluctuations have their origin in quantum 
fluctuations of the value at of the field during inflation, which translate into 
fluctuations in the number density of axions after the coherent oscillations 
of the axion held start. 

If the PQ SSB happens after inflation, i.e. if f a < Hj/2n , the present 
Universe contains many different domains that were causally disconnected 
when the SSB happened and thus have randomly distributed values of 9i. 
In this case a network of cosmic strings is generated at the moment of SSB 
of the global Upq(1) via the Kibble mechanism. This mechanism consists 
of a correlation of the held values after the SSB of a global U(l) symme¬ 
try around an axis in space and along the orbit of degenerate minima in 
internal space. This prevents SSB along the axis itself, with the conse¬ 
quent formation of a topological defect along the axis, i.e. a string (see 
e.g. [165]). These strings may become connected with walls after the ex¬ 
plicit PQ symmetric breaking, and then the whole network of strings (and 
possibly walls) decays into axions (see e.g. [166]). Understanding the evo¬ 
lution of this network is very complicated and different calculations of the 
axion density they produce differ by up to a factor of a 1000. The red hor¬ 
izontal bands in Fig. 18 show where axions constitute the whole of the DM 
according to some of these calculations (see [164] for details) for m a ~ 100 
/zeV, f a — 10 11 GeV. These masses (see Fig. 17) are not within the reach 
of current searches with cavities in ADMX (1.9 to 10 /zeV) but could be 
tested by ADMX-HF, with smaller cavities [167]. 

Finally, let us mention that axions are also relics of the pre-BBN era 
and its characteristics could change in non-standard cosmologies [168]. 


7. Sterile neutrinos as DM candidates 

The SM has three “active neutrinos” i/ 0 , i.e. neutrinos with weak interac¬ 
tions, but others without weak interactions (called “sterile” by Bruno Pon- 
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Figure 19. Prediction for sterile neutrinos which constitute the whole of the DM (red 
lines) produced non-resonantly [170] (L = 0 - the pink band is due to uncertainties in 
this production mechanism [171]) or resonantly [172] (for large lepton asymmetries L = 
0.003, 0.001 and 0.1). Also shown are limits from phase-space considerations in different 
structures (Tremaine-Gunn and Fornax Core) and non-observation of X-rays generated 
in their decay. “Pulsar-Kicks” indicates a particular region of interest (see e.g [45]). The 
red symbol corresponds to the sterile neutrino which could produce a recent potential 
line signal reported at 3.5 keV. Fig. from [176]. See [176] for explanations and references. 


tecorvo) v s can be easily added (one or more, of any mass). These are mixed 
to the active neutrinos through their common mass matrix. Considering 
only one active and one sterile neutrinos, the mass eigenstates v\ and V 2 are 
admixtures of both: i/ a = cos 6 + sin0 vi and v s = — sin0 v\ + cos 9 z^. 

Here v a ,s are interaction eigenstates and are mass eigenstates, with 
mi << m 2 = m s . The mostly sterile neutrino mass eigenstate is a good 
candidate to be WDM if its mass is m s =0(keV) [169]. This neutrino can 
be created via active-sterile oscillations, without [170] or with [172] a large 
Lepton Asymmetry L (much larger than the baryon asymmetry in the Uni¬ 
verse), and respectively be WDM or “Cool DM” (i.e. a cooler form of 
WDM, with a spectrum peaked at smaller momenta). These two types 
of sterile neutrinos are called “Dodelson-Widrow (DW)” [170] (or “non- 
resonantly produced”) and “resonantly produced”, respectively (because 
a large L produces an L-driven resonant conversion of active into sterile 
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neutrinos). DW sterile neutrinos have an almost thermal spectrum, with 
average momentum over temperature ( p/T) ~ 3.15, but resonantly pro¬ 
duced ones can easily have (p/T) ~ 1 to 2. 

A combination of lower bounds on DW sterile neutrino masses, coming 
from early structures in the Universe (called Ly-a clouds) and sub-halo 
counts in simulations of galaxy formation, leads to m s > 8 keV. Upper 
limits on the same mass coming from the non-observation of X-rays due to 
the decay v s —> 1/7 in galaxies and clusters of galaxies lead to m s < 3 keV. 
These together rule out DW sterile neutrinos as the major component of 
the DM (see e.g. [173] and [174] and references therein). The same limits 
are less restrictive when applied to resonantly produced sterile neutrinos, 
which remain valid candidates to account for the whole of the DM. Also, 
sterile neutrinos can be produced in other ways besides oscillations (e.g. 
in the decay of new scalar fields or heavier sterile neutrinos) which yield 
sterile neutrinos that can also constitute the bulk of the DM (see e.g. [45] 
and references therein). 

If they exist, sterile neutrinos are remnants of the pre-BBN era and the 
aforementioned constraints assume a standard cosmology. Their relic abun¬ 
dance could be very different in non-standard cosmologies (see e.g. [175]). 

Recently a potential weak line signal was reported at 3.5 keV [176, 177] 
which could be due to the two-body decay u s —> of a sterile neutrino with 
m s = 7 keV and sin 2 2 6 ~ 10 -10 (indicated by the red spot in Fig. 19). The 
line was found by one group in stacked observations of 73 galaxy clusters and 
in the Perseus cluster [176], and by another group in the Andromeda galaxy, 
Perseus cluster and the Milky Way center [177]. This could be a resonantly 
produced sterile neutrino [178], but other explanations in terms of atomic 
line emissions or backgrounds cannot be excluded. Further observations 
and analyses are necessary to confirm or reject this potential DM signal. 


8. Outlook 

There is no compelling observational or experimental evidence in favor of 
any of our DM candidates, such as WIMPs, axions, sterile neutrinos, pri¬ 
mordial black holes or any other of the many that have been proposed. 
It is only through experiments and observations that we will elucidate the 
nature of the DM. The next decades will be a very exciting time for DM 
research. Direct detection experiments will move to the ton-scale, indirect 
detection efforts will continue developing very rapidly, the LHC with its 
full capabilities (and possibly the next particle collider) will reframe what 
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we know about physics beyond the SM and thus DM candidates. 

The importance of the possible payoff of these searches is enormous. A 
confirmed detection of a DM candidate would open the doors to the age of 
precision DM studies to determine its properties and to DM astronomy. 
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